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ABSTRACT 

We present new planetary nebulae (PNe) positions, radial velocities, and magnitudes for 
6 early-type galaxies obtained with the Planetary Nebulae Spectrograph (PN.S), along with 
derived two-dimensional velocity and velocity dispersion fields, and the a parameters (i.e. 
the number of PNe per unit luminosity). We also present new deep absorption-line long-slit 
kinematics for 3 galaxies in the sample, obtained with the VLT/FORS2 spectrograph. 

We extend this study to include an additional 10 early-type galaxies with PNe radial 
velocity measurements available from the literature, including previous PN.S studies, in order 
to obtain a broader description of the outer-halo kinematics in early-type galaxies. These data 
extend the information derived from stellar absorption-line kinematics to typically several and 
up to ~ 8 effective radii. 

The combination of photometry, absorption-line and PNe kinematics shows: i) a good 
agreement between the PNe number density distribution and the stellar surface brightness in 
the region where the two data sets overlap; ii) a good agreement between PNe and absorption- 
line kinematics; iii) that the mean rms velocity profiles fall into two groups, with part of 
the galaxies characterized by slowly decreasing profiles and the remainder having steeply 
falling profiles; iv) a larger variety of velocity dispersion radial profiles; v) that twists and 
misalignments in the velocity fields are more frequent at large radii, including some fast ro- 
tator galaxies; vi) that outer haloes are characterised by more complex radial profiles of the 
specific angular momentum-related Ar parameter than observed within 1 R^; vii) that many 
objects are more rotationally dominated at large radii than in their central parts; and viii) that 
the halo kinematics are correlated with other galaxy properties, such as total B-band and X-ray 
luminosity, isophotal shape, total stellar mass, V/cr, and a parameter, with a clear separation 
between fast and slow rotators. 

Key words: Galaxies: general - galaxies: haloes - galaxies: elliptical and lenticular, cD - 
galaxies: kinematics and dynamics 
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de los Muchachos on the island of La Palma, of the Instituto de Astroffsica 
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1 INTRODUCTION 

The dynamics of galaxies provide fundamental information on their 
origin and evolution. Our knowledge of the kinematics and dynam- 
ics of early-type galaxies is mainly based on measurements of the 
first 4 moments of the line-of-sight velocity distribution (LOSVD) 
from stellar absorption lines. These measurements are generally 
confined to within 1-2 effective radii {Re) and can be obtained 
with long slits or integral-field units . Detailed LOSVD measure- 
ments for early-type galaxies (e.g., Bender et al. 1994: Fi s hell 19971 : 
IPinknev et alj2003l : [Emsellem et al. 2004) have provided us with a 
general picture of their dynamics. The higher-order moments are 
needed to obtain good constraints on the anisotropy and mass dis- 
tribution (Gerhard 1993; Merritt 1993). Most of the studied ob- 
jects appear to be isotropic or slightly radially an isotropic sys- 
tems jKronawitter et alj200(]| : ICapDellari et al ]l200i) . Some galax- 
ies have revealed the presence of c old stellar discs whose kinemat- 
ics stand out only after l-LSR^ jRix et al.ll 19991) . 

Kinematics of early-type galaxies are related to their isopho- 
tal shape. Deviations of the isophotes from a perfect ellipse 
are parameterised by the amplitude of the cos 40 term in a 
Fourier expansion o f the isophote radius in polar co ordinates (e.g. . 
Bender etal. 1988). Observational evidence (e.g., Ben dej|l98a : 
Kormend y & D jorgovski 1989) implies that discy ellipticals (04 > 



0) have significant rotation with V/tr > 1 and may generally be 
axisymmetric; boxy ellipticals (04 < 0) exhibit no rotation, have 
a range of values of V/cr including strongly anisotropic systems 
(V/cr << 1), may be triaxial, and are in general more massive 
than discy ellipticals. These relations between isophotal s hape and 
galaxy kinematics motivated I Kormendy & Bendeil j 19961) to revise 
the Hubble classification scheme for early-type galaxies using the 
^4 parameter, which is related to intrinsic galaxy properties rather 
than the apparent ellipticity, which is related to the galaxy's orien- 
tation on the sky. 

A re lated modification to th is scheme has been recently pro- 
posed by Emsellem et alj feOO T*), taking advantage of the advent 
of integral-field units. These allow for tw o-dimensional maps o f 
the LOSVD moments typically out to IRc ( lEmsellem er^l2004h . 
Early-type galaxies have been divided into two distinct classes: 
slow and/oi/ rotators, according to the stellar angular momentum 
they possess per unit of mass. The two classes have different dy- 
namical properties: slow rotators appear to be more massive sys- 
tems, nearly round with a significant kinematic misalignment, im- 
plying a moderate degree of triaxiality, and span a moderately large 
range of anisotropies; fast rotators appear to be rather flattened 
systems, without significant kinematic misaligrm ients, nearly ax- 
isymm etric and span a larger range of anisotropies dCappellari et al.l 
l2007h . 

In parallel, the arrival of large photometric surveys such as the 
Sloan Digital Sky Survey has consolid ated findings of a b imodal 
colour distribution in local galaxies (e.g., lBaldrv et alj2006l) . differ- 
entiating between a blue cloud of mostly starforming spiral galaxies 
and a red sequence of mostly non-starforming early-type galaxies. 

Some of the main aims of galaxy formation studies are to un- 
derstand the mechanisms that allow galaxies to evolve from the 
blue cloud to the red sequence, and to differentiate between pro- 
cesses that form early-type galaxies that are fast rotators or slow 
rotators. Numerical simulations suggest that red-sequence galax- 
ies are formed by mergers of galaxies in the b lue cloud, fol- 
lowed by a quenching of the star formation (e.g., Cattaneo et al.l 
l2006l : lFaber et alJl2007l : lKang et al.ll2007l : lRomeo et al.ll2008l) . The 
less-luminous fast rotators with discy isophotes are preferentially 



formed through a series of minor mergers with less massive com- 
panions. On the other hand, the more-luminous slow rotators 
with boxy isophotes, are thought to form through a violent major 
merger between galaxies of similar mass (e.g., Naab et al, 
Naab & Burkert 2003), o r through multiple o r hierarchical mergers 
(iWeil & Hernquist,199^ : lBurkert et al]|2008l) . 

These formation mechanisms are complicated further by the 
presence or absence of gas during the merger event, which plays a 
fundamen tal role in the final kinemat i c structure of the merger rem- 
nant ( e.g., iBarnes & Hemauisl|[T99^ : lNaab etaPl 20061 : Iciotti et all 
l2007h . Mergers producing red-sequence galaxies dominated by ro- 
tation appear to be gas-rich (wet), while gas-poor (dry) mergers 
produce red-sequence galaxies dominated by random mo tions (e.g., 
iBoumaud et al.i2005l : ICox et al.ll2006l : iNaab et alj200a) . 

However, observations show that many of the most mas- 
sive early-type galaxies were already in place by z ~ 2 
dvan Dokkum et al. 2004; Treu et al. 2005), and the evolution of the 
galaxy luminosity function since redshifts of z ~ 1 argues against 
a significant contribution of recent dissipationless dry mergers to 
the for mation of the most m assive early-type galaxies in the red se- 
quence (Scarlata et al. Therefore it is clear that the processes 
forming elliptical galaxies are not yet completely understood. 

Numerical simulations of galaxy formation in a cosmological 
context predict particular radial profiles for the total and d ark mat- 
ter distributions (e.g., Dekel et al. 20ol: lNaab et alj|2007l) . for the 
Vja ratio (e.g., Abadi et al. 2006), the angular momentum, orbital 
distribution and isophotal shape (e.g., .Naab et al. ,2006i ), depend- 
ing on which cosmology or merger type is assumed. The evalua- 
tion of these quantities from observations through dynamical mod- 
els is fundamental to probe galaxy formation theories, numerical 
simulations and cosmological scenarios. Unfortunately the picture 
is complicated by the presence of unknown variables such as the 
three-dimensional shape of the galaxy and the orbital distribution. 
Therefore only detailed information on the LOSVD out to large 
radii can help to disentangle possible scenarios in galaxy forma- 
tion. 

In spiral galaxies (or generally in gas-rich systems), the distri- 
bution of dark matter can be determined by measuring the kinemat- 
ics of neutral or ionised gas, which can be easily observed at large 
distances from the centre. This allowed for instance the discovery 
of the presence of dark matter in galaxies in the 1970s. Only a few 
early-type galaxies contain large gas rings with which the mass dis- 
tribution can be traced out to several R^ (e .g. iBertola et al.ll 19931 ; 
lOosterloo etai] |2002; Weii mans et al ] |2008l) . In the large majority 
of ellipticals or lenticulars, the mass determination is more diffi- 
cult because absorption lines from a stellar spectrum can be mea- 
sured with a sufficient signal-to-noise ratio only up to 2 effective 
radii. Studies based on integrated light spectra provide evidence 
for dark matter only in a fraction of ellipticals, with the inferred 
mass profiles being nearly isothermal to the limit of the data and 
the dark matter contributing ~ 10 - 50% of the mass within R^ (e.g . 
iKronawitter et alj|200(]| ; lOerhard et ai]|200ll ; iThomas et ai]|2007l) . 
For other systems absorption-line data are not extended enough to 
provide conclusive evidence for dark matter (see, e.g., the detailed 
recent analysis of iDe Lorenzi et alj l2008al la). The apparent di- 
chotomy may be related to the fact that the logari thmic gradient of 
the mass-to-light ratio correlates with luminosity ( INaDolitano et al.l 
I2OO5), suggesting a link between the structural parameters and the 
total mass of ellipticals. 

In order to extend the kinematic information for early-type 
galaxies to larger radii, alternative kinematic tracers have been 
identified to compensate for the rapid fall-off in the stellar light at 
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radii larger than IR^, such as globular clusters and planetary nebu- 
lae (PNe). In particular, a lot of effort has been focused on PNe in 
the last decade: they are generally believed to trace the main stel- 
lar population of elliptical galaxies, and their relative bright [O III] 
emission line allows them to be easily detected at large radii from 
the centre, making them an ideal tool for kinematic studies. 

A dedicated instrument was installed at the William Her- 
schel Telescope: t he Planetary Nebulae Spectrograph (PN.S, 
iDouglas et al. I I2OO2I) . In the last seven years, our team has under- 
taken a long-term observational campaign aimed at measuring the 
kinematics of a selected sample of ellipticals and SOs, with the 
principal aim to quantify the dark matter content in those systems. 
An initial sample of round (EO-El) and bright (S < 12.5) galax- 
ies has been selected, and later enlarged to include more flattened 
objects (< E5) and SOs. A series of papers has been published, 
presenting the initial results of the PN.S survey on early-type 
galaxies and their d y namical analysis fc omanowsky et al. 2003; 
Douglas et al. '2007'; 'No ordermeer et alTTzOOSi ; iDe Lorenzi et al.l 
2008h ; .Napolitano et al...2008l) ). 

To provide a general overview of the outer kinematics (R > 
1.5-2 Re) of early-type galaxies and fundamental constraints on 
formation scenarios, we have used the literature to collect all the 
existing PNe kinematics in ellipticals and SOs, amounting to 10 
objects in total. This includes galaxies previously observed with the 
PN.S as well as galaxies observed by other authors using different 
telescopes and techniques. 

In this paper we also present new PNe measurements for 6 
galaxies obtained with the PN.S. We will refer to them as sample 
A, while we will divide the other 10 galaxies from the literature into 
sample B (6 galaxies with more than 80 PNe radial velocity mea- 
surements) and sample C (4 galaxies with less than 80 PNe radial 
velocity measurements) in the rest of the paper. The total sample 
consists of 16 galaxies, and their principal characteristics are sum- 
marised in Table[T] In Section|2]we present the PN.S. observations, 
in Section [3] we compare the PNe spatial distribution to the stel- 
lar surface brightness. In Section|4]the smoothed two-dimensional 
velocity and velocity dispersion fields are derived and compared 
to the corresponding stellar kinematic data. In Section [5] we pro- 
vide notes on the individual galaxies. In Section |6] we discuss the 
general kinematic properties of early-type galaxies at large radii 
obtained from the combination of stellar and PNe kinematics for 
all the galaxies in the sample, and we close with a summary of our 
work in Section|7] 



2 OBSERVATIONS AND DATA REDUCTION 

For galaxies in sample A, we present new PNe positions, radial 
velocities and magnitudes, obtained as part of the PN.S elliptical 
galaxy survey. For these sample A galaxies, data reduction, PNe 
identification, radial velocity and magnitude measurements are per - 
formed using the procedures described in 'Douglas et al ] j2007h . 
Catalogues with positions, velocities and magnitudes of PNe in 
galaxies of sample A are shown in Table |2] Heliocentric systemic 
velocities derived from PNe are in agreement within the errorbars 
with the values tabulated in either the RC3 or NED catalogues. In- 
strumental magnitudes are shifted in order to match the bright cut 
off of the luminosity function with the expected apparent magni- 
tude m' = M'^fjQ-j -1-25-1-5 log D, where D is the distance in Mpc 
given in Table [T] and M^„n7 = -4.48 is the absolute magnitude of 
the cut-off according to lCiardullo et alj ([l989b). 

The number of PNe detected using the PN.S is over 100 for the 



majority of the galaxies, with a remarkable number of 457 objects 
for NGC 4374. This is a major improvement when compared to the 
results of previous techniques and instruments (typically <50 de- 
tections). In the following sections, we use only PNe whose veloci- 
ties are within 3(T of the mean velocity of their neighbours, accord- 
ing to the "friendless algorithm" described in lMerrett et alj ( l2006h . 



3 PNE SPATIAL DISTRIBUTION 

Figure [T] shows the spatial distribution of the PNe in sample A 
galaxies. For the PNe positions in the galaxies of samples B and 
C, we refer to the references listed in Table[T] 

The first property of the PNe spatial distribution we would like 
to investigate is whether or not the number density of the PNe fol- 
lows the stellar surface brightness. This has already been tested for 
many galaxies in the sample (see reference list in Table [TJ, find- 
ing generally a good agreement between the surface brightness and 
PNe counts. Q 

From Figure [T] we see that to a first approximation the two- 
dimensional distribution of PNe is symmetric. NGC 821 and NGC 
4564 are possible exceptions in that they may show an excess of 
PNe at large radii on their SE side. Unfortunately, the number of 
PNe in those regions is too small to determine whether the asym- 
metry in the distribution is statistically significant. 

In the following sections we present a comparative study be- 
tween PNe density number density and stellar surface brightness 
for the galaxies in sample A and B. For the galaxies in sample C, 
the PNe samples are too small to perform this analysis. 

In Section lSTI we describe the stellar surface brightness and in 
Section lT2l we compare it to the radial PNe number density. 



3.1 Photometric data of the stellar component 

Radial profiles of the stellar surface brightness, position angle, el- 
lipticity 6, and 04 shape parameter were obtained from the combi- 
nation of different datasets in the literature. Mean values are shown 
in Table |3] together with the sources of the photometry. 

Usually we considered HST photometry in the innermost 4 - 
5", ground based observations for r > 10-11" and their average 
values in between. 

From the e and 04 radial profile we derived their characteristic 
values by averaging the values within the ran ge 2 ■ FWHM , „j„^ < 
R < 1.5 -j?^ as done in lBender et all ll988b . |Peletier et all ( 1 1 9901) 
and lHaoetaLl l l2006l) . As weights in the average procedure, we use 
the ratios F/E where F is the flux measured at that position and E 
is the error in the (e or (34) measurement. Errors in < e > and < 04 > 
are computed by measuring the scatter of the measurements in the 
same spatial region. 

The surface brightness radial profiles p were fitted with i) a 
Sersic profile and ii) a Z?'''* plus exponential disc profil e, using the 
robust Levenberg-Marquardt method implemented by iMore et al.l 



In s ome cases this could b e complicated by population effects. For ex- 
ample, ISambhus et al] j2006l) found a subpopulation of bright PNe in NGC 
4697, which are geometrically and kinematically peculiar. To discover such 
effects in the PN distribution requires both a large PN sample and accurate 
PN magnitudes, depending on how clear the signature is. Here we only con- 
sider how well the entire PN population agrees with the surface photometry. 
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Figure 1. DSS images of galaxies in sample A with the PNe positions marked as green crosses. The field-of-view is 13 X 13 arcminutes. North is up and East 
is left. 
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Table 1. Sample galaxies 



Name Type D cz Bj PA A'™? Rlast Log Lx Reference 

(Mpc) (kms-l) (") D (') Log(ergs-l) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



NGC 821* 


E6 


22.4 ± 1.8 


1735 


11.72 


39 


25 


123(4) 


6.8 


40.33 




NGC 3377* 


E5 


10.4 ± 0.4 


665 


11.07 


41 


35 


151(3) 


10.0 


< 39.60 




NGC 3608* 


E2 


21.3 ± 1.4 


1253 


11.69 


40 


75 


87(5) 


6.8 


40.01 




NGC 4374* 


El 


17.1 ± 0.9 


1060 


10.01 


53 


135 


450(7) 


6.9 


40.83 




NGC 4564* 


E6 


13.9 ± 1.1 


1142 


12.05 


22 


47 


49(1) 


7.5 


< 39.85 




NGC 5846* 


EO 


23.1 ± 2.1 


1714 


10.91 


53 


70 


123(1) 


6.0 


41.65 




NGC 1023* 


SO 


10.6 ±0.8 


637 


10.08 


46 


87 


183(20) 


10.8 


39.6 


Noordemeer et al. (2008) 


NGC 1344 


E5 


18.4 ± 2.5 


1169 


11.24 


46 


165 


194(3) 


6.7 


< 39.48 


Teodorescu et al. (2005) 


NGC 3379* 


El 


9.8 ± 0.5 


889 


10.18 


47 


70 


186(5) 


7.2 


39.54 


Doudas et al. (2007) 


NGC 4494* 


El 


15.8 ±0.8 


1344 


10.55 


53 





255(12) 


7.6 


40.10 


Napolitano et al. (2008) 


NGC 4697 


E6 


10.9 ±0.7 


1236 


10.07 


66 


70 


535(0) 


6.6 


40.12 


Mendez et al. (2001) 


NGC 5128 


SO 


4.2 ± 0.3 


547 


7.30 


255 


35 


780 


8.4 


40.10 


Peng et al. (2004) 



NGC 1316 


SO 


20.0 ± 1.6 


1793 


9.40 


109 


50 


43 


4.3 


40.87 


NGC 1399 


El 


18.5 ± 1.4 


1447 


10.44 


42 


110 


37 


8.3 


41.63 


NGC 3384 


SO 


10.8 ±0.7 


704 


10.75 


50 


53 


68 


2.4 


< 39.52 


NGC 4406 


E3 


16.0 ± 1.0 


-244 


9.74 


91 


130 


16 


4.0 


42.05 



Arnaboldi et al. (1998) 

Saglia et al. (200O) 
Tremblav et al. ( 1995) 
Arnaboldi et al. (1996) 



Notes - Col. 1 : Name of galaxy. Galaxies marked with * have been observed with the PN.S. Col. 2: Morphological type of galaxy according to the NASA/IPAC 
Extragalactic Database (NED). Col. 3: Distance of galaxy from Tonry et al. ( 2001) derived usi ng the surface brightness fluctuations method. Distance moduh 
are shifted by -0.16 mag to take into account the new Cepheid zero point of lFreedman et al ]t.2001) as done by Jensen et al. 1 2003). Col.4: Galaxies' helio- 
centric systemic velocity liste d in NED. Col. 5 : Total B magnitude corrected for extinction and redshift listed in d e Vaucouleurs et all )l99lL RC3 hereafter) . 
Col. 6: Effec tive radii given bvlBlakeslee et alj bOOlli Col. 7: Galaxy photom etric position angle as given by RC3, except for NGC 1399 iSagha et al.l2000t) . 
NGC 3379 te apaccioU et aljl990t) and NGC 4494 iNapolitano et alj2008l) . Col.8: Total number of PNe with measured radial velocities. In parentheses we 
give the number of sources excluded by the friendless algorithm (not included in the listed total PNe number), for those galaxies for which we applied our 
adaptive kernel smoothing technique (see Section l4!2l . Col. 9: Maximum distance of the PNe detections from the galaxy centre. Col. 10: Logarithmic X-ray 
bolometric lumino sity according to O'Sullivan et al. (2001), except for N GC 1023, whose X-ray l uminosity was converted from 0.5-2 keV measurements of 
iDavid et al] i2006l) to bolometric luminosity following the prescription of l O'Sullivan et al]j200 l'). Col. 11: Reference for the PNe data. The horizontal lines 
separate the galaxies with new PNe catalogues presented in this paper (sample A), from the other galaxies whose PNe data are available in the literature, 
with more than 80 PNe radial velocities (sample B) and with less than 80 PNe radial velocity measurements (sample C). 



( I1980I) . The actual computation was done using the mpfit algo- 
rithm implemented by C. Markwardt under the IDL environment 0. 
We would like to note that the aim is not to study in detail the pho- 
tometric components of the galaxies, but instead to obtain a general 
description of the mean radial profile of the stellar light. In the cases 
where HST data were not available, we excluded the innermost 1" 
or 2" from the fit, to reduce possible contamination due to seeing 
effects. Results of the fit are shown in Table[3] 

3.2 Comparison between PNe spatial distribution and 
surface brightness profiles 

In order to compare the stellar surface brightness and the PNe num- 
ber density, we have to take into account the incompleteness in the 
PNe detections, e.g. the number of undetected PNe because of their 
low signal-to-noise ratio. The noise in the field has four main con- 
tributions, the sky surface brightness (constant all over the field), 
the detector readout noise (constant all over the field), the galaxy 
background (variable over the field) and the presence of foreground 
stars (variable over the field). 

To measure how these four noise sources affect the PNe detec- 
tion and correct the PN number counts we proceed as follows: 

^ The most recent version of the mpfit algorithm can be found at 
[http://cow.physics.wisc.edu/~craigm/idl/fitting.htmll 



Table 2. PNe catalogue for galaxies in Sample A 



Name 




R.A. 


Dec 


Vo 




"15007 


PNS-EPN- 


hh 


:mm:ss 


deg:mm:ss 


km s" 


1 


mag 


(1) 




(2) 


(3) 


(4) 




(5) 


NGC82 1-001 


02: 


08:01.6 


+ 10:58:52.4 


1670 ± 


21 


28.8 


NGC82 1-002 


02: 


08:05.5 


+ 10:58:45.2 


1703 ± 


21 


28.1 


NGC82 1-003 


02: 


08:06.7 


+ 11:02:03.4 


1742 ± 


20 


27.7 


NGC82 1-004 


02: 


08:10.4 


+ 11:00:00.4 


1725 ± 


21 


28.6 


NGC82 1-005 


02: 


08:13.0 


+ 10:58:22.4 


1823 ± 


22 


29.7 


NGC82 1-006 


02: 


08:14.0 


+ 10:57:13.7 


1781 ± 


21 


28.7 


NGC82 1-007 


02: 


08:15.8 


+ 10:59:22.3 


1669 ± 


21 


28.3 


NGC82 1-008 


02: 


08:16.1 


+ 10:59:20.0 


1832 ± 


21 


28.3 


NGC82 1-009 


02: 


08:16.2 


+ 10:56:31.0 


1629 ± 


21 


28.7 


NGC821-010 


02: 


08:16.4 


+ 11:00:17.5 


1681 ± 


21 


28.8 


NGC5846-124 


15: 


06:46.9 


+01:34:32.2 


1971 ± 


22 


28.1 



Notes - The complete table is published as Supplementary Material 
in the online version of this article. PNe catalogues are also available 
at http://www.strw.leidenuniv.nl/pns Col. 1: PNe name according to the 
lAU standards (PNS-EPN-<galaxy name>"<sequential number>). 
Col.2: Right ascension (J2000). Col.3: Declination (J2000). Col.4: Ra- 
dial velocity in the heliocentric reference system. Col. 5 : Apparent m5oo7 
magnitude. Magnitudes have an error of ±0.4 mag. 
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KGC B21 NGC 3377 KGC 360B 
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R [msec] li [.rcBec] R [„csec] 



NGC 4 374 NGC 4564 NGC 5846 




100 aOO 300 400 100 200 300 400 

H [arcsec] H [arcret;] 

Figure 2. Surface density profiles of PNe and stars. Open diamonds: Surface brightness of the stellar component as function of semi-major axis, according 
to the references in Table [T] Solid lines: Extrapolation of the stellar surface brightness profile with the blue line representing the Sersic law and green line 
representing the plus exponential decomposition. Filled circles with error bars: scaled logarithmic PNe number density (see text for details). 



• We first generate artificial images, which mimic the sky and 
CCD readout noise in our observations and then we populate them 
with simulated point-like sources of different magnitudes. 

• We compute the planetary nebulae luminosity function 
(PNLF) of the simulated objects in the artificial images. In particu- 
lar, we evaluate the magnitude mgoa at which the loss of objects is 
significant. m8o% is defined in such a way that for m < mgo%, more 
than 80% of PNe are recovered. 

• Now we consider the real, observed images and we insert a 
set of artificial point-like sources in them with a distribution in 
magnitude given by the observed PNLF, but brighter than mgo%. 
In this way we are sure (at a 80% confidence level) that our arti- 
ficial sources would not be undetected because of the background 
noise (sky plus CCD readout noise). 

• We detect the artificial sources in the real images again, and 
since we know their positions, we can compute the radial complete- 
ness factor cji, i.e., the inverse of the fraction of simulated objects 



recovered at different distances from the galaxy centr^ Since we 
inserted only sources with m < mgo%, the loss of objects is at- 
tributed only to the galaxy continuum and the stellar trails contam- 
ination rather than the sky background or detector noises. 

• Finally, we group the observed PNe in elliptical annuli, ori- 
ented along the galaxy photometric major axis. We compute N^/A, 
the number of PNe brighter than mso% in each annulus (corrected 
for the radial completeness factor c> evaluated at the annuli cen- 
tres), divided by the area A of the annulus. We then compute the 
logarithmic PNe number density ppNe(R) = -2.5 log,Q (Nc/A). 

In Figure |2] we show the comparison between the major axis 
stellar surface brightness profile and PpnAR) for galaxies in sam- 



Distance from the centre is computed taking into account the mean galaxy 
ellipticity: = l(\ - e'f- + 7^, where Y is aligned along the major axis. 
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pies A and 

• pPNeiR) is shifted by an arbitrary constant (different 
for each galaxy) to match the stellar profile. 

For NGC 1344 and NGC 4697, which are not part of the PN.S 
dataset, we could not evaluate the completeness correction with 
simulated sources because the reduced images were not available. 
This affects mostly the innermost regions of NGC 4697, in which 
the inner PNe density profile in Figure [2] falls slightly below the 
stellar surface brightness profile. 

The conclusion we can derive from this comparison is that the 
stellar surface brightness and PNe counts agree well for the galaxies 
where the two sets of data overlap in radius. In the cases where the 
stellar surface brightness measurements are not extended enough 
in radius to ensure an overlap with the PNe data, the PNe number 
density still follows the extrapolation of the surface brightness fits 
at larger radii. 



4 PNE KINEMATICS 

In this section we analyse the kinematics of the sample galaxies 
using their PNe radial velocities and compare them to long-slit 
absorption-line kinematics. 

In Section|4T|we define the "folded catalogue", which we will 
use in Section l4~2l 14.31 and 14.41 to measure the global velocity and 
velocity dispersion fields of the PNe system and their associated er- 
rors. In Section |43] we extract velocity and velocity dispersion ra- 
dial profiles from the PNe data and compare them to stellar profiles, 
to test the consistency of the two sources of kinematic information. 

This analysis is done only for galaxies in samples A and B for 
which sufficient data are available. 



4.1 Folded catalogue 

To decrease the statistical noise resulting from a low detection num- 
ber, we assume that the galaxy is point-symmetric in phase space. 
Therefore each point in phase space, (x, y, v) has a mirror counter- 
part (-X, —y, -v). Here x and y give the positions of the PNe on the 
sky, and are centred on the galaxy with y aligned along the galaxy 
photometric major axis, and v gives the radial velocities of the PNe 
corrected for the galaxy's systemic velocity (therefore v = at 
X = 0, >' = 0). By the union of the original data set and the mirrored 
one, we virtually double the number of data points. This technique 
has often been adopted in the past (e.g., Arnaboldi et al. 1998; 
IPeng et alj|2004b . Hereafter, the term "folded catalogue" refers to 
the PNe catalogue obtained by the union of the original one and 
the mirrored counterpart, and we use this one to map the 2-D kine- 
matics of the PNe system. The reliability of the point symmetry 
assumption is tested a posteriori by repeating the calculations us- 
ing only the original catalogue and comparing the two results. We 
note that in almost all cases the computed quantities based on the 
original and folded PNe samples agree well; only in the case of 
NGC 4374, the rotation field shows asymmetries of the order of 25 
km s ' (Section [5T4l . 



^ NGC 5128 is not included because of the lack of an extended surface 
brightness radial profile. 



4.2 Construction of smoothed two-dimensional velocity and 
velocity dispersion fields 

The computation of the two-dimensional velocity and velocity dis- 
persion fields from the folded catalogue was done using an adap- 
tive kernel smoothing technique, which improves on the method 
described in Peng et alj 120041) . For comparison, two examples of 
unsmoothed velocity fields are shown in Appendix IbI 

As stated in Section [3] we use only PNe whose velocities are 
within 3cr of the mean velocity of their neighbours. Generally, the 
3tr outliers we rejected are few, and they make no difference to the 
output velocity and velocity dispersions fields. The only exceptions 
are the two galaxies NGC 3377 and NGC 3608 in which small vari- 
ations in the resulting two dimensional fields are seen when outliers 
are included. However, the outliers in these two cases were more 
than 5cr away from the galaxy systemic velocity and therefore we 
are confident that their exclusion is justified. 

At every position {xp, yp) on the sky we computed the velocity 
and velocity dispersion by: 



V{xp,yp) 



d-(xp,yp) 



HiVi ■ Wi.p 



(< > - < V >^ 



(1) 



2, vf • Wi.p 



-V(xp,yp)~ -AV^ 



(2) 



where Vi is the ;-th PN velocity and AV^ is the measurement error 
(~20 km s"'), which is a combination of instr umental error related 
to the PN.S (-17 km s '. lMerrett et alj|200d) and the accuracy in 
determining the planetary nebula position (a few km s"'); w,./> is 
the distance-dependent weight for the i-th PN, defined using the 
Gaussian kernel: 



= exp 



2k{xp,ypY 



(3) 



where k is the kernel amplitude and D, = -^(x,- - xp)- + (y,- - yp)^ 
is the distance of the ;-th PN from {xp,yp). 

Errors on V and a are obtained using Monte Carlo simula- 
tions, which will be discussed in Section |431 

The weighting procedure depends on the distance D from 
{xp,yp) and on the amplitude k of the kernel. The latter is a mea- 
sure of the spatial resolution at which we would like to investigate 
the kinematics. Large values of k will lead to smoother profiles of 
the velocity and velocity dispersion fields, highlighting the general 
trend but suppressing kinematic structures at small scales. Con- 
versely, smaller values of k will allow a better spatial resolution, 
with the risk of amplifying the noise pattern due to the low number 
of PNe and the measurement errors. The optimal k must therefore 
be chosen in order to find the best compromise between spatial res- 
olution and noise smoothing. It can be different at each point in the 
observed field, since it depends on the number density of points, the 
velocity gradient and the velocity dispersion we want to resolve. 

To determine k, we proceeded as follows. At each position on 
the fields, we define k to be linearly dependent on the distance Rm 
of the M - th closest PN. We arbitrarily chose M = 20 but tested 
in a range 10 < M < 60 finding no significant differences in the 
final results, even for those galaxies with a relatively low number 
of PNe. 

k{x,y) = A ■ RM(x,y) -h S = A V(x - xmP + (y - yuf + B (4) 
where xm, ym are the coordinates of the M-th closest PN to (x,y). 
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Table 3. Photometric properties of 11 sample galaxies. 



Name Band /^^ n ;jg 

NGC [mag arcsec"^] ["] Im^g arcsec"^] 

(1) (2) (3) (4) (5) (6) 



ff,.e ;io < e > < ^4 ■ 100/fl > Reference 

["] [magarcsec-2] ["] [%] 

(7) (8) (9) (10) (11) (12) 



821 


B 


23.0 + 0.1 


39.8 + 2.0 


4.7 + 0.2 


22.0 + 0.1 


20 + 2 


224 + 0.2 


29+ 1 


0.38+0.01 


0.70 + 0.1 


(1).(2) 


3377 


B 


22.6 + 0.6 


54 + 4 


5.2 + 0.3 


21.5 + 0.2 


26 + 3 


22.1+04 


40 + 3 


0.50 + 0.01 


0.9 + 0.9 


(1),(2).(3) 


3608 


B 


25.28 ± 0.09 


157 + 10 


7.0 + 0.1 


21.9 + 0.1 


19+1 


22.5 + 0.1 


53 + 3 


0.19 + 0.02 


-0.15+0.09 


(2),(3) 


4374 


V 


23.1 ±0.2 


142+16 


8.0 + 0.6 


20.62 + 0.05 


37 ± 1 


22.4 ± 0.08 


123 + 3 


0.13+0.02 


-0.4 + 0.1 


(4) 


4564 


B 


22.1 ±0.1 


33.8 ±2 


3.1+0.2 


21.42 + 0.05 


13.3+0.5 


20.33 + 0.02 


12.90 + 0.08 


0.39 + 0.11 


0.7 + 0.9 


(1) 


5846 


V 


29 ± 1 


2903 + 192 


12 + 2 


22.52 + 0.04 


68 ±2 


24.26 + 0.07 


379 + 16 


0.07 + 0.05 


0.0 + 0.1 


(5) 



1023 


B 


214 + 0.1 


60 + 2 


3.9 + 0.2 


20.6 + 0.1 


31+3 


21.2 + 0.1 


62+ 1 


0.39 + 0.05 


0.54 ±0.1 


(6) 


1344 


V 


21.5 + 0.1 


50 + 2 


4.1+0.1 


21.27 + 0.07 


42 + 2 


21.9 + 0.5 


34 + 3 


0.31 +0.03 


0.1 ±0.2 


(7) 


3379 


B 


22.1+0.1 


47.0 + 0.2 


4.7 + 0.2 


22.25 + 0.08 


57 + 1 


17.2 + 0.5 


0.1 +0.5 


0.09 ± 0.02 


0.2 ±0.1 


(8).(9) 


4494 


V 


21.63 + 0.07 


49+1 


3.3 + 0.1 


21.74 + 0.06 


49 + 2 


22.1+0.5 


30 + 7 


0.14 + 0.01 


0.2 ± 0.2 


(10) 


4697 


B 


20.60 ± 0.04 


66+ 1 


3.53 + 0.06 


20.97 + 0.04 


74+1 


20.1+0.3 


26 + 3 


0.42 ± 0.04 


1.4 + 0.2 


(2).(ll) 



Notes - Col. 1 : Galaxy name. Galaxies in sample A are listed in the upper part of the table. The lower part of the table contains the sample B galaxies except 
NGC 5128. Col. 2: Photometric band. Col. 3: Surface brightness at Rg . Col. 4: Scale radius Rs determined from the Sersic fit. Col. 5: Sersic index. Col. 6: 
Bulge brightness at R^b for the /?''''+ exponential disc fit. Col. 7: Bulge effective radius for the /?''''+ exponential disc fit. Col. 8: Central surface brightness of 
the exponential disc for the exponential disc fit. Col. 9: Scale radius of the exponential disc component fit. Col. 10: Weighted mean value for ellipticity 

in the range 2 ■ FWHM seeing < R < \.5 ■ Re- Col. 11: Weighted mean value for t he 134 parameter in t he range 2 ■ FWHM,eemg < R < 1.5 ■ Rg. Positive 
values con'espond to discy isophotes, negative values to boxy isophotes. From Emsellem et alj fcOOvt). In that p aper the error on 04 is n ot gi ven, therefore 
we ass ume d an arbitrary error of . 1 . Col. 12: Reference for the photometric data sets: (l) lGout jfroor et al.l il994l): ( 2) | Lauer et a l] i2005l ); (3)|jedrzeiewskil 
h987l): (4)lKormendy et"ai]j2008l): (5) j Kronawitter et al. (2000); (6) Noorderme er et alj <2008l) : (T'l lsikkema et alj jloOTi) : (8) ICapaccioh et al.l h990l) : (9) 
iGebhardt et al.l ilOOd) : dO) lNat?olitano et al. l ,200&) : (1 D .De Lorenzi et al., 1,2008^1 



Our definition of can be used also to define the local PNe num- 
ber density p(x, y): 



M 



(5) 



which can be used to relate k to the local PNe density p, by com- 
bining Equation |4] and Equation|5] 



k(x,y} = AJ — + B 
I np 



(6) 



The way k is defined through Equation |4] (and its equivalent 
Equation [6} allows it to be smaller in the innermost regions where 
the spatial density is bigger and larger in the outer regions, where 
the PNe density is smaller. 

In order to select the most appropriate constants A and B, we 
built simulated sets of PNe (with the number density resembling 
the observed one) with radial velocities distributed according to a 
chosen velocity gradient and velocity dispersion, which mimic the 
observations. The artificial sets were processed with the adaptive 
kernel procedure using different values of A and B until the sim- 
ulated input velocity field was recovered. Therefore, for a given 
density, velocity gradient and velocity dispersion we have the most 
reliable values for A and B to use which give the best compromise 
between kinematic resolution and smoothing (see Table |4j. 

We show in Figure [3] the smoothed two-dimensional veloc- 
ity and velocity fields for galaxies in sample A derived with this 
technique. The two-dimensional fields of the sample B galaxies are 
presented in Appendix B. 



4.3 Errors in the velocity and velocity dispersion fields 

Errors on the derived smoothed velocity and velocity dispersion 
fields are computed by means of Monte Carlo simulations. 

For each galaxy, we built 100 datasets of PNe with simulated 
radial velocities at the same positions as in the observed (folded) 



Table 4. Typical parameters and typical errors for the smoothed, two- 
dimensional velocity and velocity dispersion fields 



Name 


A 


B 


< AV > 


< Act > 






arcsec 


km s"' 


km s"' 


(1) 


(2) 


(3) 


(4) 


(5) 


NGC 0821 


0.24 


10.68 


30 


20 


NGC 3377 


0.57 


-0.08 


25 


20 


NGC 3608 


0.08 


33.56 


40 


30 


NGC 4374 


0.97 


22.72 


25 


20 


NGC 4564 


0.24 


10.68 


20 


20 


NGC 5846 


0.00 


60.00 


50 


30 



Notes - Cols. 2-3: Values of A and B used in the kernel smoothing pro- 
cedure (see Equation 13) as determined from the simulations. Cols. 4 - 
5 typical error on the two-dimensional velocity and velocity dispersion 
fields as determined from Monte Carlo simulations (see Section l43] for 
details). 

dataset, to mimic the observations. The radial velocity for each sim- 
ulated object was calculated from the observed two-dimensional 
smoothed velocity field by the addition of a random value. The ran- 
dom value must resemble the observed velocity dispersions and as- 
sociated measurement errors and therefore it was chosen randomly 
from a Gaussian distribution centred at and with dispersion equal 
to 



(7) 



where a is the velocity dispersion measured at that position and AK 
is the velocity error. 

These simulated datasets were processed with the same reduc- 
tion script and parameters as the observed ones. The statistics of the 
simulated velocity and velocity dispersion fields give us the error 
associated with any position in the observed fields. Average values 
are given in Table |4] 

Two-dimensional error fields are not shown, but errors derived 
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Figure 3. Smoothed two-dimensional velocity (left panels) and velocity dispersion (right panels) fields of galaxies in sample A from PNe data. Spatial scale 
and orientation are given in the panels. The photometric major axis as given in Table[T]is aligned along the vertical axis. Crosses represent the locations of the 
PNe, while the colours represent the values of the smoothed velocity (or velocity dispersion) field at those points. The colour scale is given at the bottom of 
each panel. The dashed ellipses are located at 2 effective radii, which are listed in Table[T] Two-dimensional fields of the sample B galaxies are presented in 
Appendix IaI 



from them will be shown in the plots of the kinematics extracted 
along the kinematic major axis (Section [4. 3. it and along the pho- 
tometric major and minor axes (Section|43]2). 



These simulations are used also to investigate two main as- 
pects: i) if the measured rotation is significant; and ii) if the mis- 
alignment between kinematic (PA kin) and photometric (PA phot) 
major axes or the twisting of the velocity fields are significant. 



4.3.1 Testing the rotation 

We extracted the velocity curve V,„ajor along the kinematic major 
axis of the PNe system from the smoothed two-dimensional veloc- 
ity field. To find its direction, we fitted the PNe velocities with a 
simple rotation model: 



Vpn((I>pn) = V„,a.v ■ cos(0™ - PAkin) 



(8) 



where (ppi^ is the position angle of the PN on the sky, the constant 
V„,oj- measures the amplitude of rotation and PA^jn is the position 
of the kinematic major axis. Angles are measured on the sky plane, 
starting from North going counterclockwise. 
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Figure 3. Continued 



Then we created 100 simulated data sets of PNe at the same 
positions as the observed ones, with velocity equal to zero plus a 
random error The random value was generated as in Section |43I 
from a Gaussian distribution with mean zero and dispersion de- 
pending on the observed velocity dispersion and the measurement 
error (see Equation|7](. 

We built the two-dimensional velocity field from the simu- 
lated catalogues using the same procedures and the same parame- 
ters adopted for the real catalogues, and we extracted the simulated 
rotation curve along the direction of the kinematic major axis. 

Since the simulated catalogues were artificially generated to 
have their PNe velocities consistent with km s"', the rotation we 
observed on the simulated two-dimensional fields is only an arti- 
fact of the noise, as specified by the limited number of data points. 



the measurement errors and the intrinsic velocity dispersion of the 
galaxy. 

By looking at the distribution of the 100 simulated rotation 
curves, we determined the region in position-velocity space where 
the rotation is consistent with km s"' at the Icr level. In Figure|4] 
we show the velocity radial profiles measured along the kinematic 
major axis of the observed two-dimensional field together with the 
lo" confidence level of zero rotation. 

The general result is that the PNe rotation we measure is real 
(with the exception of NGC 5846). The typical Icr range in which 
the velocity is consistent with zero is around 20 - 30 km s"', de- 
pending mostly on the number of detections and the velocity dis- 
persion value. 
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Figure 3. Continued 



4.3.2 Testing misalignment and twisting 

In this section we investigate how precisely we are able to identify 
the kinematic major axis of the PNe system (and thereby ascertain 
whether the misalignment with the photometric axis is consistent 
with zero or not). We also investigate the reliability of distortions 
in the velocity field, which may be interpreted as an effect of triax- 
iality. 

To do that, we measured the kinematic position angle of all 
the simulated fields computed in Section 1431 The standard devia- 
tion of these values gives us an estimate of the error in the kine- 
matic major axis determination. These errors typically range be- 
tween 10 and 25 degrees, depending on the number of PNe, and on 



the amount of rotation and velocity dispersion: errors are smaller 
in galaxies with larger number of PNe and higher Via (see also 
iNapolitano et alj|200ll) . Results for PAkin and errors are given in 
Table|5] as well as the misalignment between kinematic and photo- 
metric major axes. The kinematic and photometric major axes are 
aligned in NGC 3377 (on average) and NGC 4564, but misaligned 
in NGC 821, NGC 3608 (marginally) and NGC 4374. 

The kinematic major axis PAkh^ might also not be constant 
with radius. To test that, for each galaxy we divided the PNe dataset 
into elliptical annuli (the ellipticity of each annulus is the galaxy's 
mean ellipticity, as given in Table [3]l and calculated the position 
angle at each radial position using Equation [8] and Monte Carlo 
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Figure 4. Red line: Radial profile of observed mean velocity for the PNe system extracted along the kinematic major axis of the two-dimensional field of 
Figure[3] Red dashed lines: Icr error on the observed mean velocity. The grey area spans the region for which the rotation is consistent with (within Icr 
level), as measured from 100 simulated datasets with zero intrinsic rotation (see Section l4.3.1l for details). 



simulations to compute the errors. This is shown in Figures|5]and|6] 
for galaxies of Sample A and Sample B respectively. 

A significant twist is observed in NGC 3377. In the inner re- 
gion (R < 80") the kinematic position angle (on the receding side) 
is PAkin = 248 ± 20, while in the outer regions it decreases linearly 
to a value of PAkin = 167 ±20 ^t R = 200". In the last measured 
bin, the rotation is very low (see Section [5.1.2b and therefore errors 
in PAkin are too large to derive definitive conclusions. For other 
galaxies, the twist is consistent with within the errors. 

Together with the kinematical twisting. Figures [5] and [6] also 
show the difference between kinematic and photometric major axis, 
using PApfiOT as zero point. An interesting case is NGC 4374: the 
mean PA^in is misaligned by 86° with respect to PApmot (see dis- 
cussion in Section [5.1.4b and the photometric position angle radial 
profile changes by 100°, reaching PAkin in the halo regions. 



4.4 Results for the velocity and velocity dispersion fields 

Smoothed two-dimensional velocity and velocity dispersion fields 
are shown in Figure [3] and Appendix B for the galaxies in sam- 
ple A and sample B. They will be used in the computation of the 
Ap (proxy for angular momentum per unit mass) (Section [6.3b , the 
outer < V/cr > ratio (Section r6.2b . and in the comparison with stel- 
lar kinematics (Section r4.5.2b . 

A general result emerging from inspection of the velocity 
fields in Figure [3] and the quantitative analysis in Figure |4] is that 
rotation is present in all sample A galaxies, except in NGC 5846 in 
which the rotation is consistent with (see Section l4.3.1l ). 

In the galaxies with significant rotation, the kinematic major 
axis may or may not be aligned with the photometric major axis, 
within the error bars (see Table|5}. Both are aligned in NGC 3377 
(on average) and NGC 4564, but misaligned in NGC 821, NGC 
3608 (marginally) and NGC 4374. In the sample B galaxies, mis- 
alignment is seen in Cen A (NGC 5128) and (marginally) in NGC 
1344 and NGC 3379. Individual cases will be discussed separately 
in Section[5] 



Table 5. Kinematic position angles and misalignments. 



Name PAkin I APA I™, Twist? PA^]^"' \APA\stars 
(NGC) (deg) (deg) (deg) (deg) 

(1) (2) (3) (4) (5) (6) 



0821 149 ±25 56 ±25 NO 31 6 

3377 221 ±20 6 ±18 YES 226 11 

3608 40 ±29 35 ± 29 NO -95 10 

4374 221 ± 18 86 ±18 NO 141 6 

4564 52 ± 10 5 ± 10 NO 49 2 

5846 237 ±86 13 ±86 NO 306 56 



1023 93 ±10 6 ± 10 NO 89 2 

1344 197 ±27 32 ± 27 NO - 

3379 -71 ±40 39 ±40 NO -108 2 

4494 180 ±17 0±17 NO 

4697 254 ± 14 4 ± 14 NO 

5128 252 ±6 37 ± 6 NO 



Notes - Col. 1: Galaxy name. The horizontal line separates galaxies of Sam- 
ple A and Sample B. Col. 2: Kinematic position angle of the PNe system 
(receding side), measured from North towards East. Col. 3: Absolute differ- 
ence between PAkin and PApnoT ± 180. PApHOT is given in Table[T] In the 
case of NGC 5846, the measured rotation is consistent with km s"' ; there- 
fore the error on PAkin is very large and the measurement of PAkin itself 
is not useful. Col. 4: Twisting of the PNe kinematic major axis, if present. 
Col. 5: Kinematic position angle of the stars (receding side), measur ed from 
two dimensionai vefocity maps, if avaiiabfe iCappettari et al]|2007l) . Col.6: 
Absolute difference between PA^J^'^* and PAphot ± 180. 

Twisting of the kinematic major axis is significant only in 
NGC 3377 (FigureH). 

4.5 Comparison with stellar kinematics 

An important aspect of our analysis is to check whether the PNe 
kinematics is in agreement with the stellar absorption-line kine- 
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Figure 5. Radial dependence of the misalignment between the photo- 
metric and kinematic major axes (as derived from the PNe kinematics) 
APA = PApHOT - PAkin (±180 deg) for galaxies in sample A (black open 
diamonds). PAphqt is the constant position angle given in Table[T] Blue 
open diamonds in the NGC 3608 panel represent the kinematic position an- 
gles calculated considering only the PNe on the north side of the galaxy. 
This sub-sample is defined and discussed in Section [3713] Significant mis- 
alignment is observed in NGC 821 and NGC 4374, while a twist in the 
direction of rotation is observed in NGC 3377. For NGC 3608 and NGC 
5846 no definitive conclusions can be derived within the errors. Red dots 
show PApHOT ~ PA(R), where PA{R) is the photometric major axis radial 
profile from the referenced papers listed in Table[3] 

N10a3 N1344 N3379 N4494 N4697 N5128 
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Figure 6. Same as Figure|5]but for galaxies in Sample B. 

matics. To do that, we retrieved major and minor (where available) 
stellar kinematic data from the literature and we compared it to the 
PNe kinematics extracted along the same axis. 



4. 5. 1 Long-slit kinematics 

Stellar kinematics from long-slit or integral-field spectroscopy is 
available in the literature for the majority of the galaxies. We give 
the list of references for each galaxy in Section|5] For two galaxies 
in sample A (NGC 3377 and NGC 4374), we obtained new deep 
long-slit observations. Also the El galaxy NGC 4494, which is part 
of the PN.S galaxy sample, was observed during this run. Its PNe 
kinematics and properties togeth er with a new dynamica l analysis 
are discussed in a separate paper jNapolitano et al.|[2008h . 



Appendix |C] gives the kinematic data derived from the long- 
slit spectra. 

4.5.2 PNe radial profiles and comparison with stellar kinematics 

Radial kinematic profiles have been extracted from the PNe data in 
two different ways. 

The first method uses the interpolated smoothed two- 
dimensional fields and their errors presented in Sections 14.21 and 
14.31 The second method selects only those PNe within an angu- 
lar section aligned along a desired position angle (the same for 
which the stellar kinematics are observed). Usually the range of 
angles is between 30 and 60 degrees, depending on the number 
of PNe. Then the selected PNe are folded to positive radii (i.e., 
the receding side of the system) and grouped into radial bins con- 
taining the same number of PNe. The number of objects Nbin per 
bin ranges from 10 to 30, to reach a compromise between num- 
ber of bins and a statistically significant number of PNe in each 
bin. In each bin the weighted mean velocity Vein and velocity dis- 
persion (Tfi/iv are computed, weights are computed from the mea- 
surements errors on the PNe velocities. The derived values of Vein 
and CTpiN change by less tha t 10 km s~' if instead a ll weights are 
equal, or if the weights from lHargreaves et al.l ( 1 19941) are used. Er- 
rors are computed with the usual formulae AVbin = ctbinI ^Nbin 
and ha BIN = a bin I ^/2(Nbin - !)■ 

Both methods have the advantage that we can extract the radial 
profile along any direction we desire, and the disadvantage that it 
might be contaminated by PNe far away from the direction we are 
interested, by the use of a large kernel smoothing parameter in the 
first method or a large angular range in the second. Especially if a 
nearly edge-on disk is present, both techniques could therefore lead 
to systematically low rotation and systematically high dispersions 
on the major axis, caused by the dilution from off-axis velocities. 
However, simple simulations show that this effect is smaller than 
the error bars in the measurement, and in any case leads to an un- 
derestimate of the significance of rotation in Figure|4l 

In Figure|7]we compare the (folded) velocity and velocity dis- 
persion radial profiles of stars and PNe. In general the stellar and 
PNe kinematics agree well, with NGC 821 as the most uncertain 
case. A more detailed description is given in Section[5] 



5 NOTES ON INDIVIDUAL GALAXIES 

5.1 Sample A 
5.1.1 NGC 821 

NGC 821 is a field elliptical galaxy. It is classified as an E6 in the 
RC3 and NED catalogues, but we measured a mean ellipticity of 
0. 4. The PNe system of NG C 821 was studied for the first time 
bv lRomanowskv et al.l j2003l , 104 PNe). The declining radial pro- 
file of the velocity dispersion observed from the combination of 
stellar and PNe data was interpreted as a signature of low dark 
matter concentration in this galaxy halo, l ower than the predic- 
tion by cosmological simulations. Recently, iForestell & Gebhardll 
(2008) extended the stellar kinematic measurements to ~ 100", 
claiming a discrepancy between stellar and PNe kinematics (the 
velocity dispersions measured with PNe are lower than the ones 
measured with absorption lines). Using three-integral models they 
de rived a dark matter content that is slightly higher than determined 
bv lRomanowskv et al.l j2003h . but still lower than predicted by cos- 
mological simulations. 
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Figure 7. Comparison between long slit stellar and PNe kinematics. Each galaxy is represented by a sub-figure consisting of radial velocity (upper panels) 
and velocity dispersion profiles (lower panels) extracted along the photometric major (left panels) and minor axes (right panels). The position angles given 
are from Table [T] Open diamonds represent the stellar kineinatics while filled circles represent the PNe kinematics extracted along a cone aligned with the 
major or minor axis. The dashed lines represents the kinematics extracted along the major or minor axis from the two-dimensional field and the dot-dashed 
lines represent the related error at the l-o" level calculated from Monte Carlo simulations (see text for details). References for the absorption-line kinematics: 
NGC 821, Forestell & Gebhardt 1 2008); NGC 337 7, our VLS/FORS2 data (see Appendi xQ; NGC 3608, Halliday et al. 1 2001); NGC 4374, our VLS/FORS2 
data (see AppendixO; NGC 4564. iHallidav et alj i200lh ; NGC 5846, Kronawitter et al.l (2000, major axis) and lEmsellem et al.l <2004 minor axis). Electronic 
tables with the kinematic data plotted in this figure are available in the on-line version of this paper 



In this paper we have obtained a new PN.S catalogue with 125 
detections. To compare their spatial distribution with the stellar sur- 
face density profile, we used ground-based photometric data in the 
V-band measured by iGoudfrooii et al.l ( Il994l) . combined with the 
HST measurements (F555W, shifted artificial ly by -0.06 magni- 
tudes to match the ground-based data) given bv fLauer et alH2005l) . 
The extrapolation of the Sersic fit to the surface brightness profile 



out to X 300" agrees well within the error bars with the number 
density radial profile of the PNe (Figure^. 

The smoothed PNe two-dimensional velocity field produced 
in Section 1431 using the folded PNe sample shows a rotation of 
~ 120 km s"', with a misalignment of 56 ± 25 degrees between 
the photometric major axis {PApnor = 25, from RC3) and the PNe 
kinematic major axis (PAkin = -31, direction of the approaching 
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Figure 7. Contin ued. References for the a bsorption-lin e kinematics: NGC 1023 , IPebattista et alj i2002L major axis) and lSimien & Prugniei J 19971 minor 
axis); NGC 1344.lTeodorescu et alj i2005h; NGC 3379.IStatler & Smecker-Hand 419991); NGc"449 4 our VLS/FORS2 data (see Appendix O; NGC 4697, 
iDeionghe et alj )l99d major axis) and lBinnev et alj Jl99(l minor axis): NGC SllS. lHuTet alj J 19951 major axis), minor axis not available. 



side). The misalignment between PAphot and PAkin is confirmed 
also if we use the un-folded PNe sample. This test excludes the pos- 
sibility that the observed rotation is produced by a few background 
sources erroneously identified as PNe, which skew the velocity dis- 
tribution. 

The comparison between PNe and stellar absorption-line 
kinematics is show n in F igure [T] Long-slit data are from 
iForestell & Gebhard 3 l l2008h, which agree w ith the integral-field 
stellar kinematics from lEmsellem et^ The PNe velocities 

at ~ 70" may be systematically slightly lower than the absorption- 
line points but the two data sets still appear consistent within the 



error bars of the PNe points and the scatter of the absorption-line 
data. 

Along the photometric major axis, the stellar component sug- 
gests a plateau in the rotation curve with V — 6Q km s"' for R < 
90", while the PNe system shows a lower value of ~ 15 - 20 ± 40 
km s"' for R > 120". Along the photometric minor axis the scat- 
ter in the stellar kinematics is quite large and it is difficult to make 
a clear comparison. Nevertheless, we see an increase of the stel- 
lar rotation in the region 20" - 50" in good agreement with the 
PNe data. The galaxy thus exhibits both major axis and minor axis 
rotation and so is most likely triaxial. In particular, the outer re- 
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gions of NGC 821 {R > 50") have the positive side of rotation 
aligned with PA = 149 ± 25 (see Figure [3] and Table [5}, while 
the inner regions have have the positiv e side of rotation aligned 
with PA = 31 jCapDellarietal.ll20d7l. but see also Figure 4 in 
lEmsellem et al. I l2004l) . indicating a misalignment of 118 degrees 
between inner and outer regions. 

The fact that PNe and absorption-line kinematics show dif- 
ferent misalignments with PApuor, 56 and 6 degrees respectively 
(Table O, could be because these two quantities are computed in 
different radial intervals: the PNe misalignment is computed for 
R > 30" -40", while the absorption-line misalignment is computed 
for R < 20" (from S AURON kinematics). There is a (small) region 
of overlap between the PNe kinematics and the long-slit data in 
Figure[7] in which there is no evidence that the PNe and absorption 
line data disagree. Both sets of data may trace different regions of 
the same velocity field, if the kinematics major axis changes rapidly 
from the centre outwards. 

The combined velocity dispersion profile along both axes 
shows a strong decrease, reaching cr 50 km s"' at 250". 



5.1.2 NGC 3377 

NGC 3377 is an E5 galaxy in the Leo I group. This galaxy has 
been studied on numerous occasions. Mass models have been ob- 
tained over different radial ranges, from the very i nner regions to 
meas u re the supermassive bl ack hole mass (e.g., lOebhardt et al.l 
I2OO3I : iMagorrian et al.l 1 19981) , out to ~ I R^ to determine the 
mass-to-light ratio, o rbital structure and dynamical properties (e.g., 
Ivan der Marel & van Dokkum 2007; Cappellari et al. 200^. The 
large number of PNe obtained with the PN.S (154 detections) de- 
rived in the present paper allows us to extend the kinematic infor- 
mation out to ~ 5 Rg. 

We used ground-based photometric data in the S-band mea- 
sured by Goudfrooij et al. (1995) and [^drzejewski (1987 ) , com- 
bined with the HST measurements given by ILauer et al.l j2005l) 
(F555W, shifted artificially by -0.08 magnitudes to match the 
ground-based data). In the small overlapping region (60" < R < 
100") the PNe counts follow the stellar surface brightness within 
the error bars. In the region 100" < R < 200" the PNe counts 
follow the Sersic fit extrapolation of the stellar surface brightness. 
The last PNe data point at 300" has too large error bars to derive 
conclusions. 

The PNe system shows a twist of the kinematic major axis, 
from PAkin = 248 ± 20 at = 60" up to PAkin = 167 ± 20 at 
R = 200" . The stellar surface brightness is not extended enough to 
check whether or not this twisting is reproduced also in the pho- 
tometry. 

Our deep stellar long-slit kinematics obtained with 
FORSA'LT (see Appendix [Cj show a good agreement within 
errors with the PNe velocity and velocity dispersion extracted 
along both major and minor axes. Along the major axis, this galaxy 
shows a very steep central rotation gradient within 5" together 
with a central velocity dispersion cr = 170 km s~' most likely du e 
to the contribution of the central black hole jOebhardt et ai]|2003h . 
The stellar radial velocity reaches a maximum value V ~ 115 
km s"' within 40" and then it goes down to ~ km s"' at 210". On 
the contrary, the velocity dispersion increases from a; 60 km s"' 
at 40" to « 100 km s"' at 100" and then it goes down again to 
80 km s"'. This turnover in the velocity dispersion is observed 
both in the stars and the PNe. Therefore it appears that despite its 
classification as an E5, NGC 3377 shows disc-like kinematics in 



its inner ~ 70" (i.e., high value of the maximum rotation velocity, 
low values of the velocity dispersion). 

5.1.3 NGC 3608 

NGC 3608 is an E2 galaxy in the Leo II group. The mass mod- 
els of NGC 3608 constructed so far are confined to the inner 
regions and are based on stellar long-slit dvan der Marell Il99ll : 
IMagorrian et al 1 9981: Gebhardt et al. |2003|) or two-dimensional 



integral-field Icappellari et alj2007h kinematics. 

Radial velocities for 87 PNe were derived in the present paper. 
It has a close companion, NGC 3607, about 6 arcminutes to the 
South, whose systemic velocity (~ 935 km s"', from RC3) makes 
the two PNe systems overlap in phase space. In our PNe analysis we 
did not consider the objects closer to NGC 3607 than NGC 3608. 
This probably still leaves some contaminants in the PNe catalogue: 
in fact, there is an over-density of detections on the side of NGC 

3608 closer to its companion (see Figure[TJ; 

We used ground-based phot ometrv bv|Jedrzeiewskil il987h in 
S-band and HST photometry by iLauer et alj ( |2005|) (F555W. ar- 
tificially shifted by 0.92 magnitudes to match the ground-based 
data). Unfortunately no spatial range exists where the stellar sur- 
face brightness and PNe data overlap, so a direct comparison is not 
possible. If we extrapolate the Sersic fit of the luminous profile, we 
find a good agreement with the PNe number denisity (much better 
than the exponential disc fit extrapolation). 

The PNe velocity field is quite noisy in the outer region and 
the PNe major kinematic axis (PA = 40° ±29) is misaligned with 
the photometric one (PA = 75). Both effects are probably caused 
by the low number of detections and the possible contamination by 
the companion galaxy. 

Stellar kinemati cs (PA = 81° and PA = -9°) are taken from 
iHallidav et all l l200ll) . Again the stars and PNe do not overlap in 
radial range, so a direct comparison between the kinematics of 
two systems is not possible. Along the major axis, the stellar kine- 
matics reveal a kinem atically decoupled core (already reported by 
IHallidav et al.ll200ll) . Extrapolation of the stellar rotation at 50" 
matches the PNe data well, which show an increase of the rota- 
tion to « 100 km s"' at ~ 150" - 200". Rotation is also visible 
along the minor axis, but the scatter and error bars are too large to 
derive general conclusions. The velocity dispersion along the ma- 
jor axis decreases from the central 200 km s"' to k 100 km s"' at 
~ 40" - 50" and then shows an increase to 170 km s"' at /? = 200" . 

As pointed out before, the PNe sample of NGC 3608 might 
be contaminated by the presence of PNe belonging to NGC 3607. 
In order to quantify the contamination effect, we exclude all PNe 
on the NGC 3607 side (i.e. southern side). The remaining cata- 
logue contains therefore only 30 PNe on the northern side of NGC 
3608, to which we apply the point-symmetric reflection as in Sec- 
tion l4.1l We will refer to this folded, northern-side catalogue as the 
"N-sample". This is a conservative approach as it will eliminate 
most of the contaminants, but it will also eliminate many PNe of 
NGC 3608 and the results may suffer from low number statistics. 

In Figures [8] and |9| we present the two-dimensional velocity 
and velocity dispersion fields, and the comparison with stellar kine- 
matics extracted along the major and minor axes using the "N- 
sample". The two dimensional kinematics of the "N-sample" are 
different from those of the original PNe sample; the new velocity 
field is characterized by more regular rotation and the new velocity 
dispersion contours have a different orientation. 

The new velocity field has a kinematic position angle of 44+33 
for R < 100" and of 72 ± 29 for R > 100", with an average 
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PAxiN = 75 ± 22. Major and minor axis rotation curves are char- 
acterized by steeper gradients: they reach the maximum of the ro- 
tation VZr = 130km S-' at 70" and = 140km s"' at 60" 
before dropping rapidly to ~ -50 km s"' after 200". The radial 
velocity dispersion profile along the major axis is systematically 
lower for the "N-sample". 

Thus the outer halo kinematics of NGC 3608 depend strongly 
on the adopted PNe sample. We therefore consider the results from 
the "N-sample" as a separate case, referring to them as NGC 3608N 
for the rest of the paper. 

5.1.4 NGC 4374 

NGC 4374 (M84) is a well-studied, bright El gala xy in the Virgo 
cluster. Mass models have bee n constructed by Kr onawitter et al.l 
l l2000h : ICappellari et alJ i EoOtI) using stellar kinematics within 1 
Rf. With our deep long-slit observations and 457 PNe radial veloc- 
ity measurements obtained in this paper we are able to extend the 
kinematic information out to 6 

We used the exte nded ground-based photometry given by 
iKormendv et"ai] ilOOA which shows a very good agreement be- 
tween stellar surface brightness and PNe number density. The data 
show large variations in the position angle of the photometric ma- 
jor axis. Beyond R = 80" the photometric position angle increases 
from 120 to 180 degrees. The shift in position angle may be an 
effect of the interaction with the nearby galaxy NGC 4406. 

The two-dimensional velocity field of NGC 4374 indicates 
that the kinematic rotation axis (PAkin = 221) is almost orthog- 
onal to the photometric one (PApuoj- = 135). Simulations caiTied 
out in Section [431 clearlv indicate that the rotation along that axis 
is significant (see Figure The misalignment between kinematic 
and photometric major axes is seen also in the stellar absorption- 
line kinematics, which agree well with the PNe data. 

The mean velocity field obtained with the "unfolded" PNe cat- 
alogue (Figure[TOl( differs from the one calculated with the "folded" 
catalogue. In particular we noticed an asymmetry: the receding side 

reaches — h50 km s"', while the approaching side reaches 80 

km s"'. This asymmetry in the velocity field is confirmed also in 
the stellar kinematics, as shown in Figure [TT] 

The "unfolded" and "folded" velocity dispersion fields do not 
differ significantly from each other. They both show almost con- 
stant velocity dispersion, with small fluctuations consistent with 
measurement errors. 

The distortion of the velocity field in NGC 4374 is unlikely 
to be caused by contamination from PNe belonging to NGC 4406 
because the two galaxies have very different systemic velocities. It 
could possibly be related to the presence of a diffuse in tracluster 
light component in this part of the Virgo cluster ( Arnaboldi et al.l 
Il996l) . or to a current interaction of NGC 4374 with NGC 4406. The 
angular separation between NGC 4406 and NGC 4374 is ~ 17' ~ 
87 kpc (at a distance of 17.6 Mpc). NGC 4406 is on the approaching 
(NE) side of NGC 4374 with a systemic velocity ~ 1300 km s"' 
lower than NGC 4374 and it is closer by about 1 . 1 Mpc (see Table 

We can exclude that the observed asymmetry is caused by 
background Lya galaxies at redshift z ~ 3.1 because they are 
too few to generate such a major distortio n of the velocity field. 
In fact, according to lCiardullo et al.l ^20021) we would expect ~ 2 
Lya contaminants for our field of view (9' x 9'), redshift range 
(Az ~ 0.0087), and completeness magnitude (28 mag). 
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Figure 8. Two-dimensional velocity (upper panel) and velocity dispersion 
{lower panel) fields of NGC 3608 obtained using the "N-sample". Labels, 
scales and symbols are as in Figure|3] 



5.1.5 NGC 4564 

NGC 4564 is an E4 galaxy in the Virgo cluster. The number of de- 
tections in our data (49 PNe, with 38 outside 2 Rg) is quite low, and 
a direct comparison of their radial distribution with the stellar sur- 
face brightness (Goudfrooij et al. 1994) suffers from small number 
statistics. The radial range where both data sets overlap is small, but 
there is a good agreement within error bars with the extrapolation 
of the Sersic fit to the stellar data. 

The two-dimensional velocity field shows a steep central ve- 
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Figure 9. The same as in Figure|2] but using the NGC 3608 "N-sample". 



locity gradient parallel to the photometric major axis. The velocity 
dispersion field has a central peak of ~ 90 km s"' (~ 50") and then 

it drops down. 

^Although the major axis stellar kinematics jHallidav et al.l 

l200lh are not extended enough to have a large overlap with the 
PNe kinematics extracted along the major axis, the extrapolation 
of the stellar rotation curve agrees well with the PNe kinematics. 
The combined major axis rotation curve reaches a plateau (~ 140 
km s"' at 20") and then declines from 50" to ~ km s"' at 180". 
The velocity dispersion decreases from the central 180 km s"' to 
~ 70 km s~' at if = 20" and then declines slowly to a; 30 km s"' at 
R = 150" - 200". The NGC 4564 major axis rotation curve resem- 
bles that of an SO galaxy rather than that of an elliptical galaxy. 



5.1.6 NGC 5846 

The EO galaxy NGC 5846 is the brightest member and the cen- 
tral galaxy in the NGC 5 846 group. Ma s s mod el ling of NGC 5846 
has been performed by van der Mareil ( 1 199 ll) ; iKronawitter et al.l 
l l2000l) : ICapDellari et alj l2007h . all using stellar kinematics within 
1 effective radius. With the PNe measurements (124 detections) ob- 
tained with our PN.S instrument, we are able to extend the kine- 
matic information of NGC 5846 out to 6 according to Table [T] 
A separate paper with detailed dynamical models combining stellar 
kinematics, PNe data and X-ray observations is in preparation (Das 
et al., in preparation). 

Stellar pho to metric data in V-band is taken from 
IKronawitter et alj ( 120001) . PNe number density and stellar 
photometry agree well within the error bars. 

The two-dimensional velocity field does not show significant 
rotation. The apparent rotation of ~ 30 km s"' along PA = 57 de- 
grees seems to be consistent with within the measurement errors 
(see Figure |4). 

Long- slit stellar data along th e major ( IKronawitter et al .I2OOOI) 
and minor dEmsellem et alj|2004t) axes are also consistent with no 
rotation, even if their radial coverage is not sufficient to ensure a 
proper comparison with PNe kinematics (Figure |7j. 

The velocity dispersion along the major axis shows a slightly 
declining profile, from the central ~250 km s"' down to ~170 
km s"' at 250". 
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Figure 10. Two-dimensional velocity field of NGC 4374 obtained using the 
adaptive kernel smoothing procedure (Section [4. 2t without imposing point 
symmetry on the PNe data set. 
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Figure 11. Comparison between the stellar kinematics and the PNe kine- 
matics from the "unfolded" data set of NGC 4374. Black symbols: Stellar 
kinematics extracted along the photometric minor axis, as shown in Figure 
|2]but without folding the two galaxy sides. Dashed red line: PNe kinematics 
extracted along the same axis of the stars, using the "unfolded" PNe sam- 
ple. Red dot-dashed lines: EiTor boundaries computed with Monte Carlo 
simulations. 



5.2 Sample B 
5.2.1 NGC 1023 

NGC 1023 is an SO galaxy , the br ightest member in the eponymous 
group. Noorder meer et al studied the kinematic data of 204 

PNe in the SO galaxy NGC 1023, obtained with the PN.S. The com- 
bined velocity curve from stellar and PNe data shows a maximum 
value of ~240 km s"' at 110" followed by a remarkable decline 
to ~50 km s"' at 360" . The peak in the radial velocity of NGC 
1023 occurs approximately at the same position as the minimum in 
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the velocity dispersion. Outside 110" the galaxy is dominated by 
a more pressure-supported component, and the velocity dispersion 
increases up to 140 km s"'. 

This was interpreted as a signature of the complex evolution- 
ary history of the galaxy: the inner regions with high rotation and 
low velocity dispersion indicate a quiescent disc formation process, 
while the outer parts with declining rotation and higher velocity dis- 
persion indicate heating from a merger or strong interaction with 
the companion galaxy. 

Figure [2] shows a good agreement between the stellar surface 
brightness and PNe number density within the relatively large er- 
rors. The PNe densities are differen t from the original values pre- 
sented in iNoordermeer et alj ( l2008h because the completeness cor- 
rection was not performed there. 

The two-dimensional PNe velocity and velocity dispersion 
fields are derived using the PNe catalogue of INoordermeer et al.l 
( l2008l) and indicate rapid rotation in the central « 100" followed 
by a decline towards larger radii (SectionlAt. 

Approximate agreement is found between absorption-line and 
PNe kinematics (Figure|7]l within the error bars. The small discrep- 
ancy in the rotation curve may be caused by the high inclination of 
the galaxy disk (~ 72° accord ing to J^oorderm eer et al. 2008). Stel- 
lar kine matics are taken from Debattista et alj ^20021) for the major 
axis an di irmen &PrugmellfT997h for the minor axis. 



preted as evidence for a lower content of dark matter than predicted 
by standard cosmologi cal simulations. 

In a recent paper, lOe Lorenzi et al. IlioOSa) presented a new 
approach to the analysis of NGC 33 79 using an N-particle model 
(NMAGIC, iDe Lorenzi et alj l2007h , combining stellar and PNe 
kinematics. They generalised previous works by relaxing the as- 
sumption of spherical symmetry and demonstrated that, due to the 
strong degeneracy between radial anisotropy, intrinsic shape of the 
galaxy and mass distribution, the data do allow moderately massive 
dark matter haloes. 

The stellar surface brightness and PNe number density 
show good agreem ent (see Figure [2] confirming the analysis of 
iDouglas et al.l2007h . 

In this paper we u se the 191 PNe radial velocities published 
bv lDouglas et alj j2007h . The smoothed two-dimensional fields are 
shown in Appendix |A] Our smoothing technique is different from 
that used previously, but the general features of the kinematics de- 
rived from the two methods are consistent (in particular, the mis- 
alignment found between the photometric and kinematic major 
axes). 

Absorption-l ine kinematics along the major and minor axes 
dStatler & Smecker-H ane 1999) agree well with the PNe data, a s 
alreadv shown in lRomanowskv et al. U2003h : lDouglas et alj ( 1200% . 
The combined velocity dispersion profile decreases to ~ 60 km s~' 
at - 210". 



5.2.2 NGC 1344 (NGC 1340) 

NGC 1344 is an E5 elliptical galaxy in the Fornax cluster. We use 
the stellar long-slit kine matics and the radial ve locities of the 195 
PNe in NGC 1344 from lTeodorescu et alj ( l2005h . who estimated a 
dark matter halo of 3.8 • 10" Mq within 160" (3.5 R,) by combining 
the two datasets. The PNe number density profile agrees well with 
the stellar surface brightness i S ikke ma et al. 2007); see Figure|2] 

The PNe velocity field (see Appendix [At reaches ~ 40 km s"' 
and shows a misalignment of 32 ± 27 degrees with respect to the 
photometric major axis. The absorption-line kinematics taken from 
iTeodorescu et alj j2005l) extend to 60" along the major axis. The 
rotation velocity reaches a maximum of ~ 100 km s"' at 60". The 
PNe velocities, which are measured o utside this region show e ssen- 
tially no rotation (already reported in lTeodorescu et al]|2005h . The 
spatial overlap between the two systems is too small to ascertain 
the cause of this steep decrease in rotation. Along the minor axis 
both components show no rotation. The velocity dispersion profile 
declines from a central value of ~ 180 km s"' to ~ 90 km s"' at 
200". 



5.2.3 NGC 3379 

NGC 3379 (M105) is an El galaxy in the Leo I group, close in pro- 
jection (7') to NGC 3384. The PNe system of NGC 3379 has been 
widely studied i n the literature. Rad ial velocities are presented by 
ICiard ullo et al. h993l 29 PNe), iRomanowskv et al. (20 oi, a pre- 
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the final sam- 



liminary sample of 109 P Ne from the PN.S), _ 
( l2006l 54 PNe) and by iDouglas et all ( 12007 , 
pie of 191 PNe from the PN.S). Also the dynamical struc- 
ture of NGC 3379 has been studied by several authors, using 
models derived from stellar kinematics ( Kronawitter et al. 20001 ; 
ICappellari et alj 2007) or from the combination of stellar and 
PNe data teornanowskv et alJl2003l ; [Samurovi£ & Danzigedl2005l ; 
IDouglas et al. I l2007l) . The PNe kinematics at large radii show a 
steep decline in the velocity dispersion, which was previously inter- 



5.2.4 NGC 

NGC 4494 is an El galaxy in the NGC 4565 group. Together 
with NGC 821 and NGC 3379, it is the third galaxy studied in 
IRomanowskv et al] j2003h in which the remarkable decline of the 
velocity dispersion radial profile was interpreted in terms of a low 
content of dark matter. A new detailed analysis is performed by 
iNapolitano et alj j2008l) . using a catalogue with almost three times 
the number of det ections obtained with t h e new PN.S pipel i ne ver - 
sion described in IDouglas et"al] ilOOH) . INapolitano et alj j2008l) 
also present deep long-slit data (see Figure ICTt and new deep pho- 
tometric data. 

The two-dimensional fields presented in Appendix |A| are 
obtained using a diffe rent technique compared to that used in 
INapolitano et alj | |2008|) . The main difference between the two ve- 
locity fields is the presence of a kinematic substructure rotating 
~ 30° off the major axis. Nevertheless, the small amplitude of this 
rotation (~ 25 km s"') is consistent with zero, if we take into ac- 
count the measurement errors and the errors determined by Monte 
Carlo simulations (~ 30 km s"'), as described in Section |43] 

Absorption-line and PNe kinematics show a good agreement, 
except at ~ 60" along the minor axis where the velocity dispersion 
measure d from the PNe is highe r than that measured from the stars 
(see also lNapolitano et alj2008l) . 



5.2.5 NGC 4697 

NGC 4697 is an E6 galax y in the direction of Virgo. This galaxy 
has 535 PNe detections bv lMendez et"ai] i200lll2008l) . They con- 
structed a spherical mass model using isotropic velocity disper- 
sions and did not re quire a dark matter halo. In this data set, 
ISambhus et"^ ( |2006|) observed the presence of two distinct sub- 
populations of PNe with different luminosity functions, spatial dis- 
tributions and radial velocities. The second popu lation was partic- 
ularly prominent in the brightest PNe. Therefore iDe Lorenzi et alj 
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( l2008bl) used only a sub-sample of 351 PNe out of the Mendez 
catalogue for their new set of dynamical models, which also incor- 
porates kinematic constraints from long-slit data. They found that a 
wide range of halo mass distributions was consistent with the kine- 
matic data for the galaxy, including some with a dark matter content 
in agreement with cosmological merger simulations. 

For this galaxy we could not correct the PNe data for the ra- 
dial incompleteness because of the lack of original images (see 
Figure O. This may explain why the PNe number density profile 
falls below the stellar s urface brightness jGoudfrooii etal .11 19941 : 
iDe Lorenzi et al ]|2008bl) in the inner regions. 

We present the velocity and velocity dispersion fields of 
NGC 4697 in (App endix [At using the sub-sample defined by 
lOe Lorenzi et alj ( [2008b.) . A clear rotation is visible along the pho- 
tometric major axis, with a peak of ~ 90 km s"' around ~ \ and 

a declining profile thereafter. 

Absorption-lin e kinematics alon g the major jPeionghe et al.l 

1 1996) and minor ( BinnevetalJll990h axes show good agreement 
with the PNe kinematics. 



kinematic m easurements from their ste llar long-sli t data, and glob- 
ular c luster tossler-Patig et alj [l998h and PNe jAmaboldi et all 
Il994h radial velocities out to 100". They found a total mass of 
1.2 - 2.5 • 10'" M0, consistent wit h the results obtained from X-ray 
determinations jlkebe et al.|[l993). A more detailed a nalysis of the 
PN sample was presented in lNapolitano et al] ^20021) . They found 
that the PN dispersion profile, depending on the definition of out- 
liers, is consistent with staying approximately flat at 200km s"' but 
also with increasing to 400 km s"' outside ~ 'iRg as the globular 
cluster profile. 

The small number of PNe (37 detections) is too low to de- 
rive reliable two-dimensional smoothed fields with our procedure. 
Moreover, there is no overlap between stellar and PNe data; there- 
fore a comparison between the PNe distribution and stellar surface 
brightness or between PNe and stellar long-slit kinematics is not 
possible. 

The combined velocity dispersion decline s from ~ 370 km s~ ' 
at the centre to ~ 200 km s"' at 170" (~ 5 Re) jSaglia et al.ll2000h . 



5.2.6 NGC 5128 

NGC 5128 (Centaurus A) is a nearby merger remnant classified 
SOpec galaxy, in the NGC 5128 group. NGC 5128 is the early- 
type galax y with the larges t number of PNe detections. It was first 
studied bv lHui et al.l ( Il995h, using radial v elocities of 433 PNe out 
to ~22 kpc, and then byl Peng et al] ( l2004l ), using radial velocities 
of 780 PNe out to ~ 80 kpc. They found a total mass enclosed 
within 80 kpc of 5 - 6 • 10" Mq, depending on the particular model 
assumed for the dark matter halo. 

The more recent analysis of IWoodlev et al.l j2007b included 
the radial velocities of 320 globular clusters. They obtained a mass 
of 1.0 • 10'^ Mo (using only PNe within 90 kpc) and 1.3 • 10'- M^ 
(using only globular clusters). 

In Section |A] we show the two-dimensional velocity and ve- 
locity dispersion fiel ds of NGC 5128 , which are very si milar to the 
ones presented by Peng et al. (2004). The long-slit data lHui et al.l 
( Il995i) covers only the central region of the galaxy. 



5.3 Sample C 
5.3.1 NGC 1316 

NGC 1316 (Fornax A ) is a bright SO galaxy in the Fornax cluster. 
lAmaboldietaLHl99 8h measured the positions and radial velocities 
of 43 PNe in this galaxy, computing an enclosed mass of 2.9 -10" 
Mq within 16 kpc. Kinematic data (PNe and stellar long-slit) used 
in our work are taken from their paper. The relatively low number 
of PNe and their sparse distribution did not allow us to construct a 
reliable two-dimensional velocity field with our adaptive Gaussian 
kernel smoothing procedure. Also the comparison with the stellar 
surface brightness suffers from low number sta tistics. In addition, 
the stellar long-slit data dAmaboldi et alj[l998) does not overlap 
sufficiently in radius with the PNe kinematics to make a compar- 
ison. The combined velocity dispersion radial profile is relatively 
flat for R> Re (see Sectionl6!4t. 



5.3.3 NGC 3384 (NGC 3371 ) 

NGC 3384 is an SO galaxy in the Leo I group, close in projection 
(7') to NGC 3379. The two galaxies have comparable systemic ve- 
locities, which makes the two PNe systems' radial velocities over- 
lap. 

Mass models of this galaxy are mostly based on stel- 
lar ki nematics and are confined to the innermost regions 
(e.g., iGebhardt et al.l l2003l : ICappellari et"ai] |2007) . There are 
no large catalogues of PNe radial velocities for this galaxy: 
iTremblay et al.' jl995l. 6 3 PNe). lDouglas et alj ( |2007|. 23 PNe) and 



5.3.2 NGC 1399 

NGC 1399 is an EI galaxy in the cent re of Fornax cluste r. This 
galaxy has been modelled in detail by ISaglia et al.l ( 120001) usmg 



Sluis & Williams! ( I2006L 50 PNe). For the an alysis in this paper 
we use the catalogue "of lTremblav et alj ( Il995l) . A cross-check be- 
tw een the Dougl a s et al . an d Sluis & Williams cat alogues is done 
in iDouglas et al.l ( l2007t) . In iTremblav et al.l ( Il995l) PNe positions 
are not published, therefore a cross-check with that catalogue is not 
possible. 

The relatively low number of PNe and their sparse distribution 
did not allow us to construct a reliable two-dimensional velocity 
field or a reliab le number density profile. The major axis stellar 
long-slit data ( F ishe j[T997l) did not significantly overlap in radius 
with the PNe kinematics. This did not allow a comparison between 
the two data sets, but the combined velocity dispersion profiles is 
shown in Figure [73] 



5.3.4 NGC 4406 

NGC 4406 (M86) is an E3 galaxy in the Virgo cluster, near in 
projection to NGC 4374. Surface-brightness fluctuations (Table[Tll 
place this galaxy at a distance of I . I ± 1 .4 Mpc closer than NGC 
4374. The X-ray emitting hot gas envelopes of both galaxies ap- 
pear to be interacting (E. Churazov, private communication). NGC 
4406 has a radial velocity of about -250 km s"' and is falling to- 
wards M87 from behi nd. The catalogue of PNe radial velocities 
dAmaboldi et alj[l996l) consists only of 16 detections, but this has 
been sufficient to constrain the velocity dispersion profile out to 1 1 
kpc(~ l.5Re). 

The relatively low number of PNe and their sparse distribution 
did not allow us to construct a reliable two-dimensional velocity 
field. Also the comparison with the stellar surface brightness suffers 
from low number statistics. 

The stellar long- slit data teenderetai]! 19941) do not overlap 



Kinematic properties of early-type galaxy haloes using PNe. 



21 



sufficiently with the PNe kinematics for a comparison. The com- 
bined velocity dispersion ranges between 200 and 250 km s"' in 
the central ~ Q.SRg, and then declines to 96 km s"' at 142". 



6 SAMPLE PROPERTIES : RESULTS AND DISCUSSION 

In this section we analyse the general properties of the whole 
galaxy sample (Table [TJ, combining the information from stars 
and PNe. In particular, we look for correlations between kinematic 
properties (such as the V/cr ratio and the velocity dispersion radial 
profile), photometric properties (such as the total luminosity, UV 
emission or the shape of the isophotes) and the ratio of PNe num- 
ber to luminosity {a parameter). All quantities determined for this 
purpose in the following sections will be presented in Table |6] 



6.1 The a parameter 

One important aspect of the comparison between PNe and the stel- 
lar su rface brightness is the evaluation of the a parameter jjacobvl 
1 1980 ). which specifies the number of PNe associated with the 
amount of light emitted by the stellar population. 

For the galaxies in the PN.S database we compute more 
specifically ffs.i.o, the ratio between the number of PNe down to 
one magnitude fainter than the PNLF bright cut-off, m*, and the to- 
tal stellar luminosity in the fi-band, Lg, expressed in solar units, as 
follows: 

• We determine the number A'™^ of PNe brighter than rago* 
accounting for completeness correction factors Cr at different radii 
from the galaxy centre, mgo* and were defined in Section |T2l 

• We integrate the fitted Sersic stellar surface brightness profile 
in the 6-band over the radial range in which the PNe have been 
detected, giving us LeiRpNe)- For NGC 4374 and NGC 5846 where 
the V-band profiles are available instead, we convert the result into 
the B-band using the total extinction-corrected colour (B - Vo)t as 
reported in the RC3 catalogue: (S- Voh = 0.94 for NGC 4374 and 
(B - Voh = 0.96 for NGC 5846. 

• We compute aBji,'-mm as the ratio Nj'^^ I LgiRpNe)- Uncer- 
tainties in the radial extent of the PNe detection range, in the Sersic 
fit extrapolation and in the completeness corrections give us an er- 
ror estimate for a. 

• Finally, we scale the measured ci'B,,„«-mgo5!, to Q's.i.o by integrat- 
ing over the PNLF. We have 



10.0 



ttfiJ O - 0!B,n 



r" F(m' ,m')dm' 

Jin* 

J,„. ^ (™ > )am 



(9) 



where F(m',m') is the analytic expression of the PNLF 
( ICiardulloetal.lll989bh . which depends only on the cut-off mag- 
nitude m* observed in our galaxies. 

For galaxie s not incl uded in our PN.S database we use the 
aB,i.o as listed bv lBuzzoni et al. (2006). 

It is known that a is related to the UV emission in galax- 
ies (IBuzzoni et al ] l2006h . We therefore compare the as.i.o values 
to the UV excess of the sample galaxies. The latter is measured 
using the total FUV (1344-1786A) magnitude from the GALEX 
database (where available) and the extinction-corrected total V 
magnitude from the RC3 catalogue. GALEX FUV magnitudes are 
corre cted for extinction u sing the relation Af^[/j/ = 8.376 • E(B - V) 
as in IWvder et al. I j2005l) . Values for the colour excess E{B - V) 
are taken from the NED database. The dgj o parameters and the 
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Figure 12. Correlation between the ag i o parameter and FUV - V colour, 
measured from total extinction-corrected magnitudes FUV (from GALEX) 
and V (from RC3). Filled circles: Sample galaxies for which a is calculated 
from our PN.S data, according to the prescription given in Se ctionl6. 1 1 Oven 
diamonds: Sample galaxies for which a values are taken from lBuzzoni et al] 
I i2006ii . 



FUV - V colours are listed in Table[6] In Figure [T2lwe plot ffg i o 
against the GALEX FUV - V colours. The figure shows a clear 
anti-correlation: galaxies with a larger UV excess, which are the 
massive ellipticals, also have a smaller number of PNe per unit lu- 
minosity, qtbi 0. 

Using stellar population models, iBuzzoni et al. I ( l2006l) inter- 
preted this trend as a consequence of the mean post-AGB (PAGB) 
core mass being smaller in massive ellipticals, which is a result of 
a higher rate of mass loss. If the mean PAGB core mass falls below 
Mere ^ 0.52 Mq one may expect a larger fraction of Horizontal- 
Branch (HB) stars to follow the AGB-manque channel, i.e., the 
stars move directly onto the high-temperature white dwarf cool- 
ing sequence after leaving the HB, thus m issing the AGB and 
PN phases entirely dGreggio & Renzinill 19901) . A larger fraction of 
AGB-manque stars is consistent both with a lower value of o-gj o 
and a strongly enhanced galaxy UV emission, as is observed for 
the more massive elliptical galaxies in the current PNe sample. 



6.2 The average V/cr ratio for stellar and PNe systems 

The average value of the ratio V/cr in a galaxy has traditionally been 
used to describe the relative importance of rotation and anisotropy 
for dynamical equilibrium and the shape. Usually, it has been cal- 
culated from long-slit kinematics using the formula 



V/cr = V,„aJcro 



(10) 



where Vi„„y is the maximum value of the observed radial velocity 
and (To is the average value of the velocity dispersion within 0.5 
(e.g. Bender et al. 1994). 

iBinnevI ( i2005h revisited the calculation of the average V/cr, 
adapting it for two-dimensional data. The new formulation is in 
terms of the sky-averaged values of V- and cr', weighted by the 
surface de nsity. For SAUR ON integral-field data, this was approx- 
imated bv lCappellari et al.l (12007): 



Of) 

(0-2) 



(11) 



where the sum extends over the data points in the two-dimensional 
field, on isophotes up to I Rg, depending on the radial extent of the 
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data. Vi and ct; are the velocity and velocity dispersion of the data 
point and F; is the flux associated with it. 

"Classical" measurements V,„nv/iTo (Equation llOt can be 
resca led to the new "two-d imensional" values by multiplying by 
0.57 jCappellari et alj2007h . 

• Stellar Via 

Stellar (Vyo')j,^„ a re taken from the SAURON survey (from 
ICappellari et"al]|2007h . where available. For the remaining sample 
galaxies in Table[T] we compute V„,a Jao from long-slit data using 
Equation [To] and rescale to the two-dimensional values by multi- 
plying with 0.57. Errors are taken to be 0.03 for SAURON values, 
or are computed using error propagation (on Equation [10} and are 
also rescaled by 0.57. 

• PNe Via- 

We compute V/cr for the PNe system using the smoothed two- 
dimensional fields (calculated in Section |4~2l (. summing over the 
positions of the detected PNe: 



Table 6. Measured parameters for the sample galaxies 



Name 0-^^^^ d-j„ a-5.0 V„„.,xASr V,.„.,_i.o Vn,u,5.o ^vZ"'o ^I^^PNe V/ostars m',mso% Q-s.l.o (FUV-V) MJLb M, 

(NGC) (kms"') (kms"') (kms"') (kms"') (km s"') (km s"') (mag, mag) (10"''Lb.o) (mag) (lO'" Mq) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12)' (13) (14) (15) (16) 
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56 ± 18 
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00 
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All _i_ A AO 
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11. ±4. 
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7 1+'' 
'■ -2.0 


14.7 ± 4.0 


1 m'3 


1 19 ± 29 


110 


79 


0.54 


123 ± 29 


209 
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0.59 


1 1 _i_ A r\i 
1 . 1 ± U.U3 


A 1 /I _L A AT 

U.34 ± U.U3 


T7 1 

2j.o, Z /.-5 


2.6 ± 1 .5 
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1.00 
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4.4j ± ' 


7.6± 0.1 






i3yy 


198 ± 47 


240 


198 


0.83 


213 ± 47 


243 


oil 
213 


0.88 


A t^A _i_ A 

U.04 ± U.z3 


A AA _L A Al 
U.UO ± U.U3 
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48. 1 ± 13.7 
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64 ±21 
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0.37 ± 0.03 
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3379 


46 ±21 


153 


69 


0.30 


46 ±21 
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77 


0.29 


0.18 ±0.02 


0.14 ±0.03 


25.5, 26.5 


24 ± 10 




g_4+2.8 


13.9 ±4.3 


3384 


72 ±20 


63 


72 


0.93 


146 ± 20 
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0.80 


1.20 ±0.20 


0.44 ± 0.03 




9.5 ± 2.5<*' 




■^■^-0.6 


3.7 ±0.7 


3608 


169 ± 43 


130 


169 


0.99 


191 ±43 


141 
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0.98 


0.45 ± 0.03 


0.05 ± 0.03 


27.1,28.0 


2.1 ± 1.1 


7.6± 0.1 




3608N 


107 ± 33 


102 


107 


0.55 


134 ± 33 


134 


134 


1.00 


0.65 ± 0.06 


0.05 ± 0.03 


27.2, 28.0 


1.7 ±0.8 


7.6± 0.1 






4374 


207 ± 33 
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0.96 


208 ± 33 
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0.16 ±0.01 
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26.9, 28.0 


5.0 ± 1.0 
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219 ±26 


219 
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1.00 


0.80 ±0.17 


0.09 ± 0.04 




3.2 ± 0.5<''> 


7.2± 0.1 


8.1!!^ 


52.4 ± 14.9 
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46 ± 15 
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0.44 
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66 
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0.18 ±0.02 


0.25 ± 0.02 


26.0, 27.0 


10. ± 3. 
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56 ± 15 


70 


58 


0.79 


58 ± 15 
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93 


0.36 


1.02 ±0.08 


0.58 ± 0.03 


26.3, 27.8 


1.6 ± 1.5 


6.9± 0.1 






4697 


92 ±20 


154 


92 


0.60 


106 ± 20 


180 
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0.59 


0.32 ± 0.02 


0.41 ±0.14 




3.8±0.3<°' 


7.9± 0.1 


6 7+2-^ 
"■'-1.7 


15.0 ±4.5 


5128 


69 ±69 
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89 


0.52 


118 ±25 


147 


133 


0.80 


0.48 ± 0.02 


0.13±0.11 




12.5 ± 0.6'°' 








5846 


170 ± 30 


210 


170 


0.81 


174 ± 30 


219 


174 


0.79 


0.12 ±0.01 


0.03 ± 0.03 


27.3, 28.3 


1.01 ±0.17 


6.7± 0.1 


11 1+3' 
-3.2 


52 ± 16 



Notes - Col.l: Galaxy name. NGC 3608N refer to the results obtained using half of the PNe sample, as explained in Section l5.1.3l 
Col. 2: Outermost value of the velocity dispersion as given by the fit. 
Col. 3: Velocity dispersion at \Rf, as given by the fit. 
Col. 4: Velocity dispersion at 5Re, as given by the fit. 

Col. 5: Minimum value of the velocity dispersion as given by the fit, scaled by cr, (from Col. 3). 
Col. 6: Outermost value of V,„,s as given by the fit. 
Col. 7: V„„, at XR^, as given by the fit. 
Col. 8: V„„, at 5Re, as given by the fit. 

Col. 9: Minimum value of Vrms as given by the fit, scaled by cTj ^ (from Col. 7). 

Col. 10: Average value of V/cr for the PNe component. For the data from two-dimensional field. Equation ! Ill was used. For the data from major axis kinematics Equation [To] was used and values have been 

multiplied by 0.57 to rescale them to the two-dimensional case, as done bv lCappellari et alj i2007l) . 

Col. 1 1 : Average value of V/o" for the stellar component, calculated with Eguation llOl on long-slit data and corrected by the factor 0.57, or using the value reported bv lCappellari et al H2007ll if available. 
Col. 12: Bright cut-oft' magnitude and nigo'/o computed for the galaxies in Figure|2] 

Col. 13: The 1 parameter computed using the PN.S data. Values marked by '°' are from lBuzzoni et alj {200^), from lCiarduUo et alj h989j) and from ljacobv et alj \\99^ . Total V band photometry for 
NGC 4374 and NGC 5846 has been converted into B band using {B-V) = 0.94 for NGC 4374 and (B-V) = 0.96 for NGC 5846 (values taken from RC3). 

Col. 14: Dift'erence between the extinction-corrected magnitudes in the fa r UV (from GALE X) and V-band (from RC3). The symbol - means that the galaxy is not present in the GALEX database. 
Col. 15: Mass-to-light ratio from single stellar population models listed in lNapolitano et all (,2P05). The symbol - means that the measurement is not available. 
Col. 16: Total stellar mass computed from the mass-to-light ratio (Col. 15), Bj and distance from Table[T] 
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(V/o-)p 



(12) 



This implicitly incorporates the weighting by the local stellar sur- 
face density, which was shown earlier to be proportional to the PNe 
number density if the completeness correction factor cr (interpo- 
lated for Xp,yp) is taken into account: Equation [T2] weights every 
region according to the completeness-corrected number of PNe in 
the region (see Section [T2t . Note that Equation [12] is an average 
over the whole PNe field, and therefore probes a much larger re- 
gion of the galaxy than the value for the stars, which was confined 
to within 1 Rg. 

For a few galaxies the number of PNe detections is too small for 
a reliable determination of the two-dimensional velocity fields. For 
these galaxies we use the PNe mean velocities and velocity disper- 
sions in bins along the major axis as given in the respective paper 
from which the data were taken. We then take the maximum of 
these mean velocities and the average of all the velocity dispersion 
values to obtain a major-axis V/cr value analogous to the classi- 
cal formula given in Equation [TO] but not confined to within 0.5 
ifj. These values are then rescaled to the two-dimensional case by 
multiplying by 0.57. 

In Figure[T3lwe compare these stellar and PNe V/cr, providing 
information on how the dynamics of a galaxy change from the inner 
to the outer region. From this figure we can see that many galaxies 
(nearly 50%) have {Vlcr)pNe > VIcr, while for the rest, the two 
values are similar. Figure [T3] also shows the dilferences between 
the distributions of flattened (>E4) ellipticals (red), round (<E3) 
ellipticals (black) and SOs (green). This shows that: 

(i) Thi-ee of the galaxies (NGC 1023, NGC 3384 and NGC 
4564) have large inner V/cr and even larger (V/(r)pNe in the halo. 
The first two are SO galaxies while the third exhibits SO-like prop- 
erties (see Section [5.1.5b . In these systems the (V/cr) rises from 
the centre to the disk-dominated region, like in the disky elliptical 
galaxies studied by Rix et al. ( 1999). 

(ii) For flattened (NGC 821, NGC 1344, NGC 3377 and NGC 
4697) and rotating, round galaxies (NGC 3379 and NGC 4494) in 
the sample, the inner and outer values for V/cr are equal to within 
the errors. 

(iii) Galaxies with a small V/cr in the inner parts may have 
either small (NGC 4374 and NGC 5846) or larger (best exemplified 
by NGC 5128 and with larger uncertainities, NGC 1399 and NGC 
3608) V/cr in the outer parts. 



6.3 The Ar radial profile for stellar and PNe systems 

In [Emsellem et al. I ( l2007l) , a kinematic classification scheme for 
galaxies is proposed based on the Ap profile, which measures the 
importance of rotation as a function of radius and is related to the 
angular momentum per unit of mass within R : 



Ar 



Y.R,F,^- 



(13) 



where the sum includes data points within R. Based on the Ar pro- 
file within ~ IRg, gal axies can be divided into two main groups 
( lEmseUem et alj2007l) : 

(i) Galaxies with Ar > 0.1 are defined Si&fast rotators. These 
have a small misalignment between the photometric and kinematic 
major axes, and have rising Ar profiles. 




Figure 13. Comparison between the two-dimensional averaged values of 
the VI tr ratio for the stellar and PNe systems. The 1:1 relation is shown by 
a dashed line. Black labels refer to round ellipticals (< E3), red labels refer 
to flattened (> E4) ellipticals, green labels refer to SO galaxies. 



(ii) Galaxies with Ar < 0.1 are defined as slow rotators. 
They exhibit a range of misalignments between the photometric 
and kinematic major axes and have flat or decreasing Ar profiles. 

To extend their analysis to the outer haloes of elliptical galax- 
ies (beyond several R^), we calculate the Ar profiles using the PNe 
kinematics for the sample galaxies. For galaxies in which the two 
dimensional field is available, we use Equation[T3] to sum over the 
velocity and velocity dispersion field values at the positions of the 
PNe, using the factor 1 / cr as the weight instead of the flux values 
F,. This incorporates the weighting by the stellar surface density 
as discussed in the case of Equation [T2] For galaxies in which the 
kinematics are available only along the major axis, we use Equa- 
tion[T3] where the sum is meant to be extended only along the major 
axis, Fi are extracted from the extrapolation of the stellar surface 
brightness and the 0.57 correction is applied. In Figure[74lwe show 
the profiles derived using the stellar kinematics (which we used to 
separate between fast and slow rotators, to be consistent with previ- 
ous works) and the PNe kinematics. To simplify matters, Figure[T4l 
is divided into two panels separating galaxies for which the stellar 
Ar is available from the SAURON data, from the remaining galax- 
ies in our sample. Several results emerge from this figure, reflecting 
the changes in both the rotation and velocity dispersion profiles: 

(i) In some galaxies, there is a marked change in the be- 
haviour of Ar at radii larger than 1 or 2 effective radii compared 
to that in the central regions. In the case of NGC 1023, NGC 1316, 
NGC 3377 and NGC 4494, after the initial increase of Ar (as ob- 
served in the majority of fast rotators within iR^) the profiles drop, 
and they reach the "slow rotator" region (NGC 1316 and NGC 
4494) or arrive close t o it (NGC 1023, NGC 3377). According to 
[Krainovic et al.l ( [2008l) . the majority of fast rotators have structures 
with disk-like kinematics, which are responsible for the high val- 
ues of Ar. In these four fast rotators, the observed outer decrease in 
Ar may imply that the light associated with the disk structure fades 
towards larger radii. 

(ii) In the galaxies of the sample, which are slow rotators ac- 
cording to the SAURON classification, the Ap profiles grow slowly 
to values of 0.1-0.3, thus entering the "fast rotator" regime (NGC 
3608, NGC 5846 and partially for NGC 4374). A larger increase is 
visible in NGC 4406, but the small number of PNe for this galaxy 
makes its Ap uncertain. The outward rise of Ar for NGC 1399 is 
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6.4 The shape of the velocity dispersion and V,„ 
profiles 



radial 
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Figure 14. Panel 1: Radial Ar profiles for galaxies common to our sample 
and the SAURON sample. Panel 2: Radial profiles for galaxies in our 
sample, which are not in the SAURON sample. Red and blue solid lines: 
Ar profiles computed from the stellar kinematics for slow and fast rotators 
respectively (those in panel 1 were kindly provided by Eric Emsellem). 
Black solid lines: Ar profiles extracted from the PNe kinematics. Dashed 
black lines in panel I connect the last SAURON data point with the first 
PNe data point, to guide the eye on the plot. These lines were not plotted 
in panel 2 to avoid overcrowding. The dotted line (Ar = 0.1) separates fast 
and slow rotator regions. In both panels some radial Ar profiles are shown 
with thicker lines for clarity. 



In order to obtain a general overview of the shapes of the velocity 
dispersion profiles tr {R) for our sample galaxies, we first param- 
eterised the major axis profiles with suitable functions. There is 
no universally valid function for all the observed profiles, so we 
experimented with various ad hoc functions and selected the one 
that best matched the observations for each galaxy. The fitted ve- 
locity dispersion profiles were then scaled by the effective radius 
(values from Table [TJ and normalised to the value of the velocity 
dispersion at l.QRg. We did not use the central velocity dispersion 
value for the normalisation (i.e., cr(R = 0)) because it might be ar- 
tificially boosted by the presence of a central supermassive black 
hole or an unresolved central velocity gradient. We performed a 
similar paramerisation for the rotation profiles to compute the V,.,,,, 
velocity profiles, where V^„,^ = a- + V^. 

According to the normalized V„„i profiles shown in Figure [TS] 
galaxies fall into two main groups. The larger part of the sample 
shows a slightly decreasing profile from the centre outwards. The 
second group of galaxies (NGC 821, NGC 3377, NGC 3379, NGC 
4564, NGC 4494, NGC 4697) show strongly decreasing Kn,.s pro- 
files. NGC 1023 appears to be anomalous in that its V„„j increases 
initially before falling steeply. This reflects the strong rotation in 
the disc of the galaxy combined with a steep central drop in the 
veloc ity dispersion profile (Figure |7] see also iNoordermeer et al.l 
I2OO8I) . For NGC 3608, the apparent increase in the Kmi profile is 
doubtful because of the possible contamination of PNe from NGC 
3607 and the large errors associated with the velocity and veloc- 
ity dispersion measurements (Figure|7]l; the radial profile of NGC 
3608N (see Section l5.1.3t does not show anomalies, and belongs 
to the first group. 

A greater variety of trends is visible in the velocity dispersion 
profiles, reflecting a larger variance in the individual contributions 
of ordered and random motions to the V,.„„, e.g., due to the strong 
disc component in galaxies such as NGC 3377. There are galax- 
ies in which the velocity dispersion remains high (i.e., the drop in 
velocity dispersion is less than 50% of the normalisation value); 
those, which exhibit a big drop (i.e., the last value of dispersion is 
less than 50% of the central value); and those that show an increase 
towards large radii, as in the case of NGC 1023, NGC 3377 and 
NGC 3608. In NGC 3608, the observed increase is again doubtful. 
In the case of NGC 4406, the massive drop is based on the results 
derived from a small number of PNe velocities, which show an as- 
sociated increase in rotation velocity. 



slower; this galaxy remains below the "slow rotator" demarcation 
line. 

(iii) The remaining galaxies in the sample are "fast rotators" 
throughout the radial range probed. The Ar profile either rises or 
stays almost flat throughout. 

The large radial extent of the PNe data enables us to follow the 
rotation properties of elliptical galaxies out into the halo, providing 
new constraints on the processes involved in the formation of these 
galaxies (see Section|6j6j. 

Although the kinematics of the stars within R^ as measured by 
SAURON (and most traditional long-slit studies) reflect half of the 
total stellar luminosity and mass, the same is not true for the total 
stellar angular momentum. As discussed in AppendixiDl these stars 
may represent only ~ 10% of the galaxy's angular momentum, and 
an adequate global rotational picture of the galaxy requires obser- 
vations to ~ 5Re- 



6.5 Correlations between kinematic, photometric and 
morphological properties of the sample galaxies 

In this section we study the relations between some physical prop- 
erties of early-type galaxies and their kinematic properties at large 
radii. 

In the the left panels of Figure [T6] we show how the outer- 
most values of V,-,,,,, (VnnsMST, measured from the ad-hoc fitted 
profiles of Figure [Tsll are related to the galaxy's X-ray luminosity, 
total B-band luminosity, mean isophotal shape parameter < 04 >, 
mean (V7cr)/j^, measured from the PNe data, and a^j o parame- 
ter In the right panels of Figures [T6] we show similar plots, us- 
ing the minimum value of V„„s normalized by the value at l.ORg, 

MiN/Vrms.i.o in place of V, 

Similarly, Figure [TT] shows corresponding plots with the last 
values of the velocity dispersion (cTlast' measured from the ad-hoc 
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Figure IS. Comparison between radial profiles of V,-,,,, {top panel) and ve- 
locity dispersion (bottom panel) of the sample galaxies, obtained by com- 
bining the stellar and PNe kinematics along the major axis. Profiles have 
been scaled to the effective radius and normalised to their value at 1 .OR,, . 
Colours are chosen in order to highlight the contrast between lines and thus 
better distinguish different profiles. 



fitted profiles of Figure [Tst and cr^^^ja^^, the minimum value of 
velocity dispersion normalized by the velocity dispersion at l.QR^. 

The outermost values of V„„, and velocity dispersion charac- 
terize their typical values in the galaxy halo and therefore also the 
galaxy mass, while the ratios between their minima and their values 
at \.QRf give a measure of how much these profiles fall towards the 
outer radii. 

We note that using the values of Vrms and velocity dispersion 
at a fixed radius of 5Rg in these figures, instead of their outermost 
values, gives very similar results. Also, plotting the physical param- 
eters against the logarithmic gradient of V™, or velocity dispersion 
between 1^^ and 5Rg results in similar trends as shown in the right 
parts of Figures[T6]and[T7] For reference, the Vrmsi^Re) and (T(5Re) 
interpolated from the parametric fits are given in Table |6] 

From Figures [16] and [T7] we notice that more luminous galax- 
ies tend to have larger values of V^jsjast and cTlast- This is not 
surprising given that luminosity, V„„j and velocity dispersion are 
known to be related to the galaxy mass. What is new here is that 
we are exploring a larger radial range, probing the relation for the 
galaxy halo. We notice also that galaxies with higher Kni.!.Msr and 
(Tj^j^^j tend to have only boxy profiles (04 < 0), while galaxies with 
low values of K/j7.!,LAsr and o"last have a wider range of shapes 
(-1 < £(4 < 2). This is also a reflection of the k nown trend for mas - 
sive ellipticals to be more boxy in shape (e.g., lBenderetalJll989h . 
iNapolitano et al. (2005) have also found the 04 and other galaxy 
parameters (luminosity, stellar mass and central surface brightness 
profile) to correlate with the mass-to-light ratio gradients, showing 
that there might be a link between the galaxy structural parameters 
and the total mass of the galaxies. 



In addition, galaxies which have the highest peaks in the Ar 
parameter (NGC 1023, NGC 3384 and NGC 4564, see Figure[l4l( 
have fainter B magnitudes and higher values of mean {V/cr)piyg. 
Moreover, galaxies with higher Vrms.tAST and (Tj^^^^j have smaller 
cg i values (i.e., less PNe per unit luminosity). As discussed in 
Section [611 this is probably a consequence of massive early-type 
systems harbouring a larger proportion of stars on the Horizontal 
Branch that do not enter the PN sta ge. The correla tion between 
ffs i and (X was already explored by iBuzzoni et al.l ( I2006) . using 
the central velocity dispersion measurements. 

We see in the right panels of Figures [16] and [TT] that these re- 
lations still hold if we replace V„„.,xAsr with the normalized min- 
imum of V,;,,,, (i. e. V„„j.M/jv/K»«.i.o) and if we replace cr^^^^ with 
the normalized minimum velocity dispersion (i. e. o^„in/o"io)' ^^^n 
if the scatter is larger jf). This is a consequence of more massive 
galaxies having preferentially larger values of V„n.i,LA5T> o"last ™d 
flatter profiles of V„„j and velocity dispersion. 

In Figure [18] we plot the total stellar mass M, computed for a 
subsample of galaxies as a function of V„,„,la57, V„„s,min /yrms.i.o, 
''"last' ''"min/''"io and the ffgj.o parameter. The stellar mass is com- 
puted using the total B luminosity from Table [T] and the mass-to- 
light ratio in the B-band listed in Napolitano et al. (2005) (where 
available). As expected from Figures [T6]and[T7] M, correlates with 
the cf parameter, V,-ms.LAST, and with ct^ast' although with larger 
scatter. 

The general message we learn from Figures I16I18I is that 
galaxies with higher V,ms,LAST, and higher ct^^ast ™d flatter pro- 
files (i.e., higher V„„,,.m,a./ V,,,,, , and higher cr„,N/o-,„) tend to be 
the more luminous, more massive galaxies, are more pressure sup- 
ported at large radii (i.e., VjapNg < 1), with boxy isophotes (i.e., 
fl4 < 0) and tend to form less PNe. T hese results extend the pic- 
ture described bylBenderet^flgggf) based on stellar kinematics 
within 1 Rg and of Buzzoni et al. (2006), to larger radii. 

In Figures [16] and [T7] we also differentiate between fast rota- 
tors, fast rotators with a declining Ar profile in the halo, slow rota- 
tors and slow rotators with Ar > OA in the halo. On average, fast 
rotators fall on the left side of the plots, (i.e. they have lower values 
of Vrms.LAST, o'LAST^'^d morc declining profiles) while slow rotators 
fall on the right side of the plots (i.e. they have higher values of 
^rms.LAST^ ''"tAST^'^d flatter profiles). This is a reflection in the halo 
kinematics that fast (slow) rotators are on a verage less (more) mas - 
sive and have more discy (boxy) isophotes dEmsellem et al J 20071) . 

It is also interesting that slow rotators with Ar > 0.1 in the 
outer haloes on average are located between fast and slow rotators 
in Figures [76] and [771 as if they represent a link between the two 
classes. 

6.6 Comparison with galaxy formation models 

The results presented in Sections [4] and [6] on the outer halo kine- 
matics of our sample galaxies provide new constraints for models 
of elliptical galaxy formation. Most of the merger simulation papers 
to date compare their remnants to data within an effective radius or 
so, but there are a few predictions for the kinematics at larger radii. 

Line-of-sig ht velocity fie l ds for binary disc merger remnants 
are published in jjesseit et al.l ^20071) . The progenitor galaxies in 
these simulations include pure stellar discs as well as discs con- 
taining 10% of their mass in gas. Some 3:1 mergers with gas lead 
to remnants showing velocity fields with rapid rotation (e.g., their 

^ Normalization at d.lRe or 0.5Rg gives similar results. 
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remnant 3 IGS 19), with a peak velocity at 2-3Re, similar to the case 
of NGC 1023. Strong misalignments such as observed in NGC 
821 are more characteristic of the 1:1 merger remnants (both dry 
and gas-rich, e.g., 1 lClO/1 1S8), which also include slowly rotating 
remnants with radially increasing rotational support (e.g., 11C6). 
Their models also include one which shows a ring-like depression 
of the velocity dispersion ai R ^ 0.5Re associated with a corre- 
sponding increa se in (11S2), somewh at similar to the case of 
NGC 3379 (see [Be Lorenzi et 'S l2008bh . Its velocity dispersion 
profile decreases by a factor of ~ 2 from the centre out to ~ l.SRg. 
This model originated from a merger of two disc galaxies with 
spin axis perpendicular to the orbital plane. O verall, the projected 
kinematics of the merger remnants analysed bv lJesseit et al.l ^20071) 
show a variety of features seen also in our data. A more careful 
comparison including several diagnostics (V/cr values, /Is-profiles, 
velocity fields, dispersion profiles) would clearly be profitable. 

Mean rotatio n and velocity dispersion profiles are shown by 
iNaab et alJ J ioo g) for a 3:1 merger remnant from the same set of 
simulations, out to 4Rg. The kinematics of this object are charac- 
terised by rapid rotation and a falling dispersion profile (by a factor 
~ 1.5), leading to a major axis V/cr profile increasing to V/o" - 2 
at 2-4 R^. Only the SO galaxies and NGC 4564 in Figure [13] reach 
these values. By cont rast, the M, galax ies formed in a cosmological 
setting (described in iNaab et alj|2007h rotate more slowly and are 
characterised by V/cr < 0.5 with similar values in the halo and the 
central parts (disregarding the counter-rotating core seen in one of 
the remnants). This is more typical for the slowly rotating galaxies 
in Figures[T3land[T4l 

The kinematic s of the luminous haloes of isolated galaxies 
are investigated by lAbadi et"al] j2006l) . The outer haloes of these 
systems are predominantly made of stars accreted during previ- 
ous merger events. These haloes have near-spherical triaxial den- 
sity distributions, strongly increasing radial anisotropy, and cor- 
respondingly falling line-of-sight velocity dispersion profiles (by 
more than a factor of 2 out to 0.1 times their virial radius). Our data 
probe only partially into these haloes, but it is possible that some 
of the strongly falling dispersion profiles in Figure[T5]are related to 
such accreted haloes. On the other hand, the galaxies in our sam- 
ple with slowly falling dispersion profiles or those with twisting 
(NGC 3377) or misaligned rotation (NGC 821) may not be easily 
reco nciled with these mo dels . 

lOiiorbe et al] j2007h studied the properties of elliptical-like 
galaxies generated in self-consistent hydro-dynamical simulations. 
These galaxies resemble the characteristics of slow rotators, i.e., 
are massive spheroidal systems without extended disk components 
and very low (cold) gas content. Their velocity dispersion profiles 
(normalized at IR^) are generally flat or slightly declining, with val- 
ues of o"^in/o"j„ larger than 0.7 out to 6Re- These values are con- 
sistent with what we observe for the haloes of slow rotators (see 
Figure [TtIi. Galaxies in our sample with more steeply declining ve- 
locity dispersi on profiles are generally fast rotators, which are not 
included in the lOfiorbe et al.l i 20070 simulations. 

In summary, we expect that more detailed comparisons be- 
tween data as presented in this paper and the halo kinematics pre- 
dicted by galaxy formation models will shed new light on the 
merger formation histories of elliptical galaxies. 



7 SUMMARY 

We explored the outer-halo kinematics and properties of early-type 
galaxies using the PNe data available for a sample of 16 objects. 



Out of the sample, 6 new catalogues with positions, radial velocities 
and magnitudes of PNe are presented for the first time (sample A, 
namely, NGC 821, NGC 3377, NGC 3608, NGC 4374, NGC 4564 
and NGC 5846). 

We compared the radial distribution of PNe with the stellar 
surface brightness profiles, for all galaxies in sample A, and for 
galaxies in sample B, which have extended photometry available 
from the literature (NGC 1023, NGC 1344, NGC 3379, NGC 4494 
and NGC 4697). We found that: 

(i) The PNe number density profile follows the stellar surface 
brightness profile in all the galaxies in which we have a spatial 
overlap between the stellar light and the PNe counts. 

(ii) In the galaxies in which we do not have a spatial overlap, the 
PNe counts follow the extrapolation of the stellar surface brightness 
profile. 

(iii) For the galaxies with a far UV magnitude measured with 
GALEX, we confirmed results from previous studies that the total 
number of PNe per unit luminosity (q-b_i.o parameter) is related to 
the UV colour excess (i.e., galaxies with higher UV emission tend 
to have fewer PNe). 

We then determined smoothed two-dimensional velocity and 
velocity dispersion fields using an adaptive Gaussian kernel, for 
the galaxies in sample A and in sample B (namely NGC 1023, 
NGC 1344, NGC 3379, NGC 4494, NGC 4697 and NGC 5128). 
We compared the PNe kinematics with absorption-line kinematics 
available in the literature. For NGC 3377, NGC 4374 and NGC 
4494 new long-slit data were presented. The kinematic analysis 
showed that: 

(i) Rotation of the PNe system is observed in the majority of the 
studied galaxies. 

(ii) There is a good agreement between the absorption-line and 
PNe kinematics along the major and minor axes. 

We combined the PNe and absorption-line kinematics along 
the photometric major axis for all 16 galaxies in our sample, ex- 
tracting the information from the two-dimensional fields where 
available. This allowed us to probe their kinematics much further 
out than with the use of stellar kinematics alone (usually limited to 
the innermost 1-2 Rg). This gave the following results: 

(i) In several fast-rotator galaxies, kinematic twists and mis- 
alignments are observed at large radii, which are not seen in the 
SAURON data within R, (NGC 821, NGC 3377, and marginally, 
NGC 3379). 

(ii) The average V/cr of the stellar component (which probes the 
kinematics in the inner 0.5 - 1 Re) is equal to or lower than the 
average value computed for the PNe (which probes the kinematics 
for R > 2 Rf). This indicates that for a fraction of galaxies, the 
kinematics become increasingly supported by rotation in the outer 
parts. 

(iii) The radial profiles of the Af, parameter (related to the angu- 
lar momentum per unit mass) show a more complex radial depen- 
dence when their values in the halo region are taken into account. 
While the majority of fast rotators remain so in their haloes, the A/j 
profiles of slow rotators grow slowly to values of 0.1 - 0.3, there- 
fore requiring a slight modification to the classification scheme 
proposed by Emsellem et al. (2007) (NGC 3608, NGC 4374 and 
NGC 5846). Some fast rotators however (NGC 1316, NGC 1023, 
NGC 3377 and NGC 4494), have profiles that strongly decrease 
outwards, probably due to the presence of disc-like structures that 
dominate the kinematics within 1-3 R, and then fade. 
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(iv) The normalized V,,,,, profiles show that our sample galaxies 
fall into two main groups; the first group (NGC 1316, NGC 1344, 
NGC 1399, NGC 3384, NGC 4374, NGC 4406, NGC 5128 and 
NGC 5846) shows a slightly declining profile from the centre out- 
wards, the second group (NGC 821, NGC 3377, NGC 3379, NGC 
4494, NGC 4564 and NGC 4697) shows a strongly declining V,.„„ 
profile. An exceptional case is NGC 1023, where Vrms increases 
until 2.5i?e, before falling steeply. 

(v) The radial profiles of the velocity dispersion show a variety 
of shapes. There are nearly flat profiles, in which the velocity dis- 
persion falls only by a few percent from the central values; strongly 
declining profiles, in which we observe a drop of a factor of 2 in 
the velocity dispersion towards outer radii; and a few anomalous 
galaxies, which exhibit rising profiles, again related to the presence 
of a disc. 

(vi) More luminous galaxies (brighter total B magnitude and X- 
ray luminosity) tend to have flatter Km.! and flatter velocity disper- 
sion profiles and larger values of Vrms and dispersion measured at 
the outermost observed point. This is related to the fact that more 
massive galaxies have nearly flat dispersion profiles, and generally 
higher values for the outer velocity dispersion. Moreover, there is 
evidence that more massive and luminous galaxies have on average 
a lower number of PNe per unit luminosity (smaller a) than less 
massive galaxies. 

(vii) Slow rotators have on average flatter V„nj and flatter ve- 
locity dispersion profiles and larger values of V,-,„,5,la.s r and cr^^j,^. 
Conversely, fast rotators have on average more steeply declining 
Vrms and velocity dispersion profiles, with lower values of VrmsMsr 

(viii) Galaxies with high values of f^^A.sx '^^ nearly flat velocity 
dispersion profiles (i.e., more massive galaxies) are preferentially 
boxy in shape (04 < 0) and have smaller values of Vja. Galaxies 
with small values of cr^j^^^ or declining velocity dispersion profiles 
(i.e., less massive galaxies) have a larger range of shapes (-1 < 
04 < 2) and V/cr ratios. 



These results show that a full picture of the kinematics and 
angular momenta of elliptical galaxies requires information about 
their outer velocity fields. When comparing to models of the for- 
mation of ellipticals it is important to take into account the halo 
kinematics - the dynamical timescales in the haloes are longer and 
therefore the imprint of the formation mechanisms are preserved 
more strongly. 
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APPENDIX A: TWO-DIMENSIONAL VELOCITY AND 
VELOCITY DISPERSION FIELDS OF GALAXIES IN 
SAMPLE B 

Different authors adopted their own procedures and criteria to cre- 
ate smoothed two-dimensional fields of the PNe, depending on the 
amount of available data. Here, for a homogeneous analysis, we de- 
termine the two-dimensional velocity and velocity dispersion fields 
for some of these galaxies, adopting the procedure described in 



Table Al. Typical parameters and errors for the smoothed two-dimensional 
velocity and velocity dispersion fields of the galaxies in Sample B. 



Name 


A 


B 


< AV > 


< Act > 






arcsec 


km s"' 


fan s"' 


(1) 


(2) 


(3) 


(4) 


(5) 


NGC 1023 


0.48 


1.36 


30 


20 


NGC 1344 


0.53 


-2.50 


30 


20 


NGC 3379 


0.57 


14.92 


30 


30 


NGC 4494 


0.48 


11.36 


30 


20 


NGC 4697 


0.40 


12.80 


40 


20 


NGC 5128 


0.08 


18.56 


15 


f5 



Notes - Cots. 2-3: Values of A and B used in the kernei smoothing pro- 
cedure (see Equation|4) as determined from the simulations. Cols. 4 - 
5 typical error on the two-dimensional velocity and velocity dispersion 
fiefds as determined from Monte Cai'lo simulations (see Section l431 for 
details). 

Section |4~2l We do this only for the galaxies with a sufficient num- 
ber of PNe (> 80), namely NGC 1023, NGC 1344, NGC 3379, 
NGC 4494, NGC 4697 and NGC 5128 (i.e., sample B). Parameters 
used in the kernel smoothing procedure and typical errors obtained 
with Monte Carlo simulations are listed in Table IaTI Results are 
shown in Figure IaTI 



appendix B: UNSMOOTHED VELOCITY FIELDS 

The smoothing technique described in Section 1421 is required to 
measure the mean velocity and the velocity dispersion along the 
line of sight of the PNe system. A single PN radial velocity mea- 
surement randomly deviates from the local mean velocity by an 
amount that depends on the local velocity dispersion. One can ask 
how the unsmoothed map of PNe radial velocity measurements 
compares with the smoothed two-dimensional velocity field. 

In galaxies with high Vja, we expect that the velocity map 
of individual PNe radial velocity measurements preserves the char- 
acteristics of the mean velocity field. On the contrary, in galax- 
ies with low Vja we expect it to be chaotic and the mean rota- 
tion to be hidden in the scatter of the velocity dispersion, which 
dominates the kinematics. As examples, we show in Figure lBTI the 
two-dimensional maps of individual PNe radial velocities for NGC 
1023 (a system dominated by rotation) and NGC 4374 (a system 
dominated by velocity dispersion). 



APPENDIX C: LONG-SLIT KINEMATICS 

Long-slit observations for NGC 3377, NGC 4374 and NGC 4494 
were carried out with the Very Large Telescope (VLT) at the Eu- 
ropean Southern Observatory (ESO) in Paranal (Chile) from De- 
cember 2005 to April 2006 (Proposal 76.B-0788A) in service 
mode under dark time conditions. The Unit Telescope 1 (Antu) 
was mounted with the Focal Reducer/low dispersion Spectrograph 
(FORS2), which was equipped with the 1400V Grism and the 0'.'5 
slit. Spectra were taken along both the major and minor axes. Basic 
data reduction was performed using standard ESO-MIDA^ and 

* MIDAS is developed and maintained by the European Southern Obser- 
vatory. 
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Figure Al. PNe smoothed two-dimensional velocity and velocity dispersion fields for galaxies in sample B. Symbols, orientations and scales are as in Figure 

m 
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Figure Al. Continued 
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Figure Bl. Two-dimensional velocity fields of the individual PNe radial velocity measurements for NGC 1023 (left panel) and NGC 4374 (right panel), 
without applying the adaptive kernel smoothing. 




Figure CI. Radial velocities, velocity dispersions, /13 and /14 coefficients for NGC 3377 along the major and minor axes. 
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IRAF routines. All the spectra were bias-subtracted, flat-field 
corrected by quartz lamp and twilight exposures, cleaned of cos- 
mic rays, and wavelength calibrated. After calibration, the different 
spectra obtained for a given galaxy along the same position angle 
were co-added using the centre of the stellar-continuum radial pro- 
file as a reference. The contribution of the sky was determined at the 
two edges of the resulting frames where the galaxy light was neg- 
ligible, and then subtracted. Spectra from adjacent rows had been 
binned together in order to ensure a signal-to-noise ratio of at least 
25. The instrumental FWHM measured on the comparison spectra 
is « 140 km s"'. 

The stellar kinematics (mean velocity, velocity dispersion, hj, 
and /i4 Gauss-Hermite coefficients ) was measured by means of th e 
Penalized Pixel-Fitting method bv lCappellari & Emsellenj j2004) . 
kinematic stellar templ ates were chosen fr om the Indo-U.S. Coude 
Feed Spectral Librarv jValdes et aflliooi) and the MILES library 
dSan chez-Blazg uez et al. 20061 ). arid then convolved with a Gaus- 
sian function to match the instrumental FWHM. Spectral regions 
with known emission lines (H/3, [O III]) were masked and not in- 
cluded in the fit. 

In Figure ICTI we show the measured long-slit kinematics. Ta- 
ble|C| containing the measured absorption-line kinematics is avail- 
able as Supplementary Material in the online version of this article. 



R„, as the rule-of-thumb target radius for measuring the angular 
momentum. is defined as the radius where a point mass would 
have the same specific moment of inertia as the extended mass pro- 
file. Using the formula 



' l°°m)R'dR] 



j;m)RdR) 



1/2 



(D2) 



we find very similar radii to the Rlp calculations reported above, 
e.g., = 2AR,A-1Re, \OARe for Sersic n = 2,4,6. 



APPENDIX D: ROTATION AND ANGULAR MOMENTUM 
CONSIDERATIONS 

Since the angular momentum of a star is weighted by its galac- 
tocentric radius, the outer regions of a galaxy are particularly im- 
portant for estimating the total angular momentum. One cannot in 
general observe the three-dimensional positions and velocities nec- 
essary to calculate the full angular momentum, so we will use a 
two-dimensional projected proxy. This proxy angular momentum 
can be simplified to a one-dimensional problem by adopting the 
approximation that the rotation is aligned with the photometric ma- 
jor axis, and its azimuthal depend ence is related to the m ajor-axis 
velocity by V{^) = Vmax co&{(p) (cf. lKrainovic et aljboOSh . 

The cumulative proxy angular momentum inside a radius R is 

then: 

r" 

LP{R) = n V„UR'MR')R'-dR', (Dl) 
Jo 

where the stellar surface density is 1,(R). With the simplest assump- 
tions that the galaxy light follows the R'^* law, and that the rotation 
velocity Kmax and the stellar mass-to-light ratio are both constant 
with radius, one can calculate that 9% of the total proxy angular 
momentum resides inside R,,, and 50% inside 5.4 R,,. The latter 
"angular momentum effective radius" (= Rlp), is naturally larger if 
Knax increases with R or if the galaxy light is more extended (e.g. 
Rip =12. IR^ for Sersic index n = 6). It is smaller if V,^ax decreases 
or if the light is more compact (e.g. Rlp = lAR^ for n = 2). 

These calculations show that probing well outside R^ is es- 
sential for characterizing the angular momentum of an elliptical 
galaxy. The dimensionless cumulative an gular momentum-like pa- 
rameter Ap introduced bv iEmsellem e t al. ( 2007) can also be shown 
to have the same limitation, so that the global angular momentum 
requires measurements to > 5Re- 

An alternative approach is to consider the radius of gyration. 



' IRAF is distributed by NOAO, which is operated by AURA Inc., under 
contract with the National Science Foundation. 
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Figure CI. Continued, for NGC 4374 and NGC 4494. 



